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V ■ Abstract 

'^ ' We present a field-theoretic renormalization group analysis of Abanov and Chubukov's model of 

the spin density wave transition in two dimensional metals. We identify the independent field scale 
and coupling constant renormalizations in a local field theory, and argue that the damping constant 

i-rt . of spin density wave fiuctuations tracks the renormalization of the local couplings. The divergences 

at two- loop order overdetermine the renormalization constants, and are shown to be consistent with 
^, ' our renormalization scheme. We describe the physical consequences of our renormalization group 

equations, including the breakdown of Fermi liquid behavior near the "hot spots" on the Fermi 
K^ ' surface. In particular, we find that the dynamical critical exponent z receives corrections to its 

QQ \ mean-field value z = 2. At higher orders in the loop expansion, we find infrared singularities 

similar to those found by S.-S. Lee for the problem of a Fermi surface coupled to a gauge field. 

\/^ ■ A treatment of these singularities implies that an expansion in 1/A^, (where N is the number of 

^■"■^ ■ fermion flavors) fails for the present problem. We also discuss the renormalization of the pairing 
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vertex, and find an enhancement which scales as logarithm-squared of the energy scale. A similar 
enhancement is also found for a modulated bond order which is locally an Ising-nematic order. 



I. INTRODUCTION 

There is little doubt that the quantum transition involving the onset of spin density 
wave (SDW) order in a metal is of vital importance to the properties of a variety of cor- 
related electron metals. This is amply illustrated by some recent experimental studies. In 
the cuprates, Daou et al} argued that the Fermi surface change associated with such a 
transition was the key in understanding the physics of the strange metal. In the pnictide 
superconductors, experiments^"^ have explored the interesting coupling between the onsets 
of SDW order and superconductivity. In CeRhlns (and other '115' compounds), Knebel 
et al.^ have described the suppression of the SDW order by pressure, and the associated 
enhancement of superconductivity. 

The theory of Hertz^"^ has formed much of the basis of the study of the spin density wave 
transition in the literature. The central step of this theory is the derivation of an effective 
action for the spin density wave order parameter, after integrating out all the low energy 
excitations near the Fermi surface. A conventional renormalization group (RG) is then 
applied to this effective action, and this can be extended to high order using standard field- 
theoretic techniques^. However, it has long been clear that the full integration of the Fermi 
surface excitations is potentially dangerous, because the Fermi surface structure undergoes 
a singular renormalization from the SDW fluctuations. 

Important advances were subsequently made in the work of Abanov and Chubukov^^'^^. 
They argued that the Hertz analysis was essentially correct in spatial dimension li = 3, but 
that it broke down seriously m. d = 2. They proposed an alternative low energy field theory 
for d = 2, involving the bosonic SDW order parameter and fermions along arcs of the Fermi 
surface; the arcs are located near Fermi surface "hot spots" which are directly connected by 
SDW ordering wavevector. They also presented a RG study of this field theory, and found 
interesting renormalizations of the Fermi velocities at the arcs. 

This paper will present a re-examination of the model of Abanov and Chubukov, using 
a field-theoretic RG method. We will begin in Section II by introducing the low energy 
field theory for the SDW transition in two dimensional metals, and reviewing the Abanov- 
Chubukov argument for the breakdown of the Hertz theory. Section HI will define the 
independent renormalization constants using the structure of the local field theory, and 
determine their values using the divergences in a 1/A^ expansion (where N is the number 
of fermion flavors) to two loop order. Actually, the two-loop divergences overdetermine the 
renormalization constants, but we will find a consistent solution: this is a significant check on 
the consistency of our renormalization procedure. While our renormalizations of the Fermi 
velocities agree with those of Abanov and Chubukov, we find significant differences in the 
other renormalizations, and associated physical consequences. At two-loop order, the ratio 
of the velocities scales logarithmically to zero (as specified by Eq. (3.40)), and consequently 
we are able to compute RG-improved results for a variety of physical observables (which 
differ from previous results^*^'^^): 



• The non-Fermi liquid behavior at the hotspot is controlled by the fermion self energy 
given by Eq. (3.44). 

• Moving away from the hot spot, we find that Fermi liquid behavior is restored, but 
the quasiparticle residue and the Fermi velocity vary strongly as a function of the 
momentum (py) along the Fermi surface: these are given in Eq. (3.45). 

• The bosonic SDW spectrum does not obey dynamic scaling with z = 2, but instead 
obeys the 'super power-law' form in Eq. (3.46), and the amplitude of the spectrum 
scales as in Eq. (3.47). 

Going beyond two-loops, we also explored the consequences of a strong-coupling fixed point 
at which the velocity ratio and other couplings reach finite fixed-point values. Here the boson 
and fermion Green's functions obey the scaling forms in Eqs. (3.19-3.22), and the non- Fermi 
liquid behavior at the hotspot is specified in Eq. (3.23). Moving away from the hotspot, we 
have the Fermi liquid form in Eq. (3.24), with the Fermi velocity and quasiparticle residue 
given by Eq. (3.25). 

In Section IV, we describe the structure of the field theory at higher loop order. Similar 
to the effects pointed out recently by S.-S. Lee^^ for the problem of a Fermi surface coupled 
to a gauge field, we find that there are infrared singularities which lead to a breakdown 
in the naive counting of powers of 1/A^. However, unlike in the problem of a gauge field 
coupled to a single patch of the Fermi surface^^, we find that the higher order diagrams 
cannot be organized into an expansion in terms of the genus of a surface associated with the 
graph. Rather, diagrams that scale as increasingly higher powers of A^ are generated upon 
increasing the number of loops. 

In Section V, we consider the onset of pairing near the SDW transition, a question 
examined previously by Abanov, Chubukov, Finkel'stein, and Schmalian^^"^^. Like them, 
we find that the corrections to the d-wave pairing vertex are enhanced relative to the naive 
counting of powers of 1/A^. However, we also find an enhancement factor which scales as 
the logarithm-squared of the energy scale: this is the result in Eq. (5.6). We will discuss the 
interpretation of this log-squared term in Section V. 

In Section VI we show that a similar log-squared enhancement is present for the vertex 
of a bond order which is locally an Ising-nematic order; this order parameter is illustrated 
in Figs. 22 and 23. The unexpected similarity between this order, and the pairing vertex, 
is a consequence of emergent SU(2) pseudospin symmetries of the continuum theory of the 
SDW transition, with independent pseudospin rotations on different pairs of hot spots. One 
of the pseudospin rotations is the particle-hole transformation, and the other pseudospin 
symmetries will be described more completely in Section II. 




FIG. 1. Square lattice Brillouin zone showing the Fermi surface appropriate to the cuprates. The 
filled circles are the hot spots connected by the SDW wavevector Q = (7r,7r). The locations of the 
continuum fermion fields ip^ and ^1 is indicated. 



II. LOW ENERGY FIELD THEORY 



We will study the generic phase transition between a Fermi liquid and a SDW state in two 
spatial dimensions, and our discussion also easily generalizes to charge density wave order. 
The wavevector of the density wave order is Q, and we assume that there exist points on 
the Fermi surface connected by Q] these points are known as hot spots. We assume further 
that the Fermi velocities at a pair of hot spots connected by Q are not parallel to each other; 
this avoids the case of 'nested Fermi surfaces', which we will not treat here. 

A particular realization of the above situation is provided by the case of SDW ordering 
on the square lattice at wavevector Q = (tt, tt). We also take a Fermi surface appropriate for 
the cuprates, generated by a tight-binding model with first and second neighbor hopping. 
We will restrict all our subsequent discussion to this case for simplicity. 

At wavevector Q = (vr, vr) the SDW ordering is coUinear, and so is described by a three 
component real field </>", a = x,y,z. There are n = 4 pairs of hot spots, as shown in 
Fig. 1. We introduce fermion fields (■?/'io-,'02o-)5 ^ = l...n, a =ti for each pair of hot spots. 
Lattice rotations map the pairs of hot spots into each other, acting cyclically on the index 
i. Moreover, the two hot spots within each pair are related by a reflection across a lattice 
diagonal. It will be useful to promote each field ip to have A^-flavors with an eye to performing 
a 1/A^ expansion. (Note that in Ref. 14, the total number of hot spots 2nN is denoted as 
A^.) The flavor index is suppressed in all the expressions. The low energy effective theory is 
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FIG. 2. Configuration of the i = 1 pair of liot spots, with the momenta of the fermion fields 
measured from the common hot spot at A; = 0, indicated by the filled circle. The Fermi velocities 
■iTi^2 of the ^^1,2 fermions are indicated. 



given by the Lagrangian, 

L = ^(dr^) + y (V0)' + — <? + —i4>'f 

+ Ar (CC'^L + <C'^L') (2.1) 

The first hne in Eq. (2.1) is the usual 0(3) model for the SDW order parameter, the second 
line is the fermion kinetic energy and the third line is the interaction between the SDW 
order parameter and the fermions at the hot spots. Here, we have linearized the fermion 
dispersion near the hot spots and v^ are the corresponding Fermi velocities. It is convenient 
to choose coordinate axes along directions x = 4=(1, 1) and y = 4=(— 1, 1), so that 

vt^ = {v^,Vy) , if=^ = {-v^,Vyy, (2.2) 

these Fermi velocities are indicated in Fig. 2. The other Fermi velocities are related by 
rotations, v^ = {Rn/2Y~^v^^^- 

We choose the coefficient A of the fermion-SDW interaction to be of (9(1) in A^. As a 
result, the coefficients in the first line of Eq. (2.1) are all scaled by A^ as this factor will 
automatically appear upon integrating out the high-momentum/frequency modes of the 
fermion fields. 

Before proceeding with the analysis of the theory (2.1), let us note its symmetries. Besides 
the microscopic translation, point-group, spin-rotation and time-reversal symmetries, the 
low energy theory possesses a set of four emergent SU{2) pseudospin symmetries associated 
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FIG. 3. Modification of the Fermi surfaces in Fig. 2 by SDW order with (0) ^ 0. The full lines are 
the Fermi surfaces, and the white, light shaded, and dark shaded regions denote momenta where 0, 
1, and 2 of the bands are occupied. The upper and lower lines are boundaries of hole and electron 
pockets respectively. 

with particle-hole transformations. Let us introduce a four- component spinor, 

(2.3) 




We will denote the particle- hole indices in the four-component spinor by a,/3. The spinor 
(2.3) satisfies the hermiticity condition, 

^-'{[~o]^' = ^f (2.4) 

Then, the fermion part of the Lagrangian (2.1) can be rewritten as, 

L^ = l^'M - ^4 ■ V)v&^i + l^'id. - ^4 ■ ^)^i + ^A0 • {¥,'r¥, + ¥,'r¥,) (2.5) 

Now the Lagrangian (2.5) and the condition (2.4) are manifestly invariant under, 

SU{2)e: ¥,^Ue¥, (2.6) 

with Ui - SU{2) matrices. We note that the diagonal subgroup of (2.6) is associated with 
independent conservation of the fermion number at each hot spot pair. The symmetry (2.6) 
is a consequence of linearization of the fermion spectrum near the hot spots and is broken 
by higher order terms in the dispersion. The diagonal subgroup noted above is preserved 
by higher order terms in the dispersion, but is broken by four-fermi interactions, which map 
fermion pairs from opposite hot spots into each other. Both symmetry breaking effects are 
irrelevant in the scaling limit discussed below. 

The pseudospin symmetry (2.6) constrains the form of the fermion Green's function to 
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be, 

which iinphes, 



-(*L^r,l') = S'-5,,5^,5^^,Glix - x') (2.7) 



G^(x - a;') = -G^ix' - x) (2.8) 

The corresponding expression in momentum space, Gl{k) = —G^{—k), imphes that the 
location of hot spots in the Brillouin zone is not renormahzed by the spin wave fluctuations 
in the low energy theory. 

Another important manifestation of the particle-hole symmetry is the equality of any 
Feynman graphs, which are related by a reversal of a fermion loop direction. 



A. The Hertz action 

The Hertz action is derived by working in the metallic phase, and integrating out the 
fermions in Eq. (2.1), leaving an effective theory for (p alone. In particular, the one-loop 
self-energy of the field (p is evaluated in Appendix A 1, and is given by 

U\uj,q) = ll\io = 0,q = 0) + N^\u\ + ..., ^ = - (2.9) 

ZTTVxVy 

The presence of the non-analytic term \uj\ is due to the fact that the density of particle-hole 
pairs with momentum Q and energy u scales as u. As usual, the constant piece 11° (g = 0) 
is eliminated by tuning the coefficient r. The ellipses in Eq. (2.9) denote terms analytic in 
u and q, starting with u^ and (f. These terms formally disappear when we take the cut-off 
of the effective theory (2.1) to infinity. Thus, the quadratic part of the effective action for 
the field reads 

52 = f I 9l^<^'^(-^' -^) (71^1 + ^y + k' + r^ r{k, uj) (2.10) 

At sufficiently low energies, the analytic term cu^ in the boson self-energy coming from the 
bare action, Eq. (2.1), can be neglected compared to the dynamically generated l^l term. 
Thus, at low energies the propagation of collective spin excitations becomes diffusive, due 
to the damping by the fermions at the hot spots. 

Hertz^ proceeds by neglecting all the quartic and higher order self-interactions of the field 
0, which are generated when the fermions are eliminated. This is justified if such interactions 
are local, as one can then absorb them into operators, which are polynomial in the order 
parameter and its derivatives (the simplest of which is just the operator (0^)^). The theory 



then reduces to, 

Sh = yI ^'^"(-^' -^) (^1^1 + ^' + '^"^^' "^^ + T^ / ^^^'^('^^')' (2.11) 

The quadratic part of the action (2.11) is invariant under scahng with the dynamical critical 
exponent z = 2, 

k ^ sk, cj — )■ s^w, 0(x, r) — )• s0(sx, s^r) (2.12) 

Thus the theory is effectively d + z = 4 dimensional and the quartic coupling u is marginal 
by power-counting in rf = 2. 

At one loop order, the flow of u follows easily from the conventional momentum shell 

dU 11 r, , , 

where s = e~^ is the renormalization scale. Thus u is marginally irrelevant, and flows to the 
Gaussian fixed point with m = in the infrared. This stability of the Gaussian fixed point 
has formed the basis of much of the subsequent work^'^'^^ on the Hertz theory. 



B. Breakdown of the Hertz theory 

The analysis in Section II A is valid only under the assumption that the fermion- induced 
quartic and higher order couplings of the field (p can be neglected. In fact, as observed in 
Refs. 11 and 14, this assumption is not justified in spatial dimension d = 2. Indeed, as 
shown in Ref. 14, the fermion-induced four-point vertex is given by, 

Tr"'"''iquq2,q3,q4) = AV"^'^^'^^"^(gi, ^2, gs, ^4) + permutations of 2, 3, 4 (2.14) 

J [qi,q2,q3,q4) - y~^ TT \ ^ -.£ /^ , -NX/-/ — ] — ^ — -^ /- , -NX 

^ 27rVxVy{t{uJ2 + UJ3) - Vi ■ [q2 + q3))[t[uJi + UJ2) - v^ ■ [qi + ^2)) 

(2.15) 

We see that the vertex (2.14) is highly non-local. Moreover, under the z = 2 scaling 

(2.12), we can neglect the frequency dependence in the denominators of Eq. (2.15), obtaining 

r4 ~ \^\/<f ~ C^(l); which produces a marginal interaction. Similarly, one can show that 

all the higher order fermion-induced vertices behave as T2n ~ li^l/l^^""^ ~ l^^"^"") which 

is again marginal under (2.12) when combined with the scaling of the field-strength. Thus, 

the Hertz-Millis theory has an infinite number of non-local marginal perturbations and the 

standard action (2.11) is incomplete. 
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C. RG interpretation 

An RG interpretation of the results of Section II B follows by performing a scaling analysis 
directly on the spin-fermion model (2.1). As before, we will scale the boson fields according 
to Eq. (2.12). Correspondingly, it is natural to scale the fermion momenta towards the hot 
spots, 

iji^{x,T)^s^l^ij[^{sx,s^T) (2.16) 

Here the field-strength rescaling has been chosen to preserve the spatial gradient terms in 
the fermion action. We now see that the boson- fermion coupling A in (2.1) is marginal under 
the field scalings in Eqs. (2.12) and (2.16); a similar analysis in c? = 3 would show that A is 
irrelevant. 

The marginality of A, and the infinite number of marginal couplings in Section II B 
indicate that all subsequent RG should be performed direction on the spin-fermion model 
(2.1). Further, with the scalings as in (2.12) and (2.16), we should not expand in powers of A, 
but rather analyze the theory at a fixed boson- fermion "Yukawa" coupling. A similar strategy 
was followed in Refs. 18 and 19 for the Ising-nematic transition in a rf-wave superconductor. 

An important consequence of the scalings (2.12) and (2.16) on (2.1) is that both the boson 
kinetic term (9t-0)^ and the fermion kinetic term ip'^drip are irrelevant. We may safely drop 
the boson kinetic energy. However, the fermion kinetic energy must be retained - otherwise, 
the theory does not possess any dynamics. We will return to this point shortly. Let us 
now rescale the fermion fields iJj = ^/v\ to eliminate the marginal coupling A. We define, 
?7 = 1/A and v = v/\ . Note that v has the unusual dimensions of [a;]^/^/[A;]. We drop the 
tildes in what follows. Then, 

, N.^-'.o Nr ^^ Nu.-.2,2 
^ = Y(V0)' + — 0' + — (0')' 

+ r {4ir:^'4a' + 4ir:^'^L') (2.17) 

As already remarked, the coupling constant r] is irrelevant. Thus, we take the limit 77 — )■ O"*" 
in all our calculations. In practice, r] gives the prescription for integrating over the poles of 
the fermion propagator. We will work with the action (2.17) for the rest of this paper. At 
criticality it is characterized by two dimensionless constants, 

a = ^ , u = - 2.18 

Vx 7 

and a dimensionful constant 7, Eq. (2.9), 

7 = ^^. (2.19) 



27rVxVy 
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FIG. 4. The boson self-energy at A^ = oo. The full lines represent the ^^1,2 fermions, and the 
dashed lines represent the boson ^°. 



Thus, in the critical regime, the theory (2.17) does not possess an expansion in any couphng 
constant. 



III. FIELD-THEORETIC RG 

We begin by discussing the general renormalization structure of (2.17). In the absence 
of a coupling constant, we will use the RPA based scaling (2.12) and (2.16) as the starting 
point of our analysis. Naively, one expects that this scaling is also obeyed by the N = oo 
limit of the theory and that corrections to it can be calculated in a systematic expansion 
in 1/A^. Indeed, the usual arguments would indicate that at A^ = oo, the boson self-energy 
is given by the RPA bubble in Fig. 4, Eq. (2.9), (see the Appendix Al for details of the 
calculation). Hence, the bosonic propagator 

(0"(x)0''(x')) = S^^Dix - x') (3.1) 

Sit N = oo takes the form, 

^(-) = 11 ^"^'? , , \ e— ^-^-^ (3.2) 

^ ^ N J (27r)3 7|a;| + g^ + r ^ ^ 

which respects the scaling (2.12). On the other hand, the fermion propagator 

-(C(^)^J^(^')) = S'"^S.A^,G',{x - x') 
at A^ = oo is given by its free value, 

duj<fk 1 



Gl{x) = I ^^, 1_^,-...+....^ (3.3) 
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Applying scaling (2.16) to this propagator indicates t] scales to zero; we will eventually take 
this limit, but need a non-zero rj for now to properly define the fermion loop integrals. 

As we will see later in Section IV, the A^ = oo limit in the present theory turns out to 
be much more subtle and is not given by the simple forms in Eqs. (3. 2), (3. 3). Moreover, 
the anomalous dimensions in this limit are not expected to be parametrically small. Never- 
theless, we can reasonably expect that the RG structure presented here remains valid, even 
though we are not able to accurately compute higher loop corrections to the renormaliza- 
tion constants. In addition, the difficulties with the 1/A^ expansion appear only at high 
loop order, which enables us to check the consistency of our approach to the order discussed 
below. 

With the above remarks in mind, we are ready to discuss the renormalization of the theory 
in Eq. (2.17). The theory contains five operators that are marginal by power counting 
at 2; = 2, and not related by symmetry. Two of these are eliminated by field-strength 
renormalizations, 

= ZyV., i^ = Z'J'A (3.4) 

As is conventional, we can fix Z^ by demanding that the coefficient of (V^)^ remains invari- 
ant. For fermion field, it is convenient to allow both velocities to fiow, and so we renormalize 
these as 

v. = z:v:, vy = zx. (3.5) 

The fermion spatial gradient terms are then not available to fix Z^, and we cannot use 
the fermion temporal gradient term because its coefficient rj scales to zero. Instead we 
demand the invariance of the boson-fermion coupling term to fix the fermion field strength 
renormalization; it is thus consistent to use a unit coefficient for this term, as we have done 
in Eq. (2.17). The quartic boson coupling renormalizes 

ZuZi^Z^ ^ 

u = — -2 — Ur. (3.6) 

It is also useful to track the renormalization of the dimensionless velocity ratio a in Eq. (2.18) 

Zv 

All the renormalization factors Z depend only on A^, Or, Ur and the ratio /i/A, where /i is 
a renormalization scale and A is a UV cutoff. 

An important point is that the damping parameter 7 appearing in the boson propagator 
does not have an independent renormalization constant. It is not a coupling in a local field 
theory, and only appears in certain correlation functions as a measure of the strength of the 
particle- hole continuum, as determined by Eq. (2.19). This implies that when we consider 
the renormalization of the boson propagator, the renormalization of the parameter 7 should 
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track the the renormahzations of the velocties Vx,y obtained from the renormahzation of the 
fermion propagator; in other words, the renormahzation of 7 is 
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z^zf"-- 



(3.8) 



This tight couphng between the boson and fermion sectors is a key feature of the theory 
(2.17), and a primary reason for strong couphng physics in d = 2. 



The theory (2.17) contains two relevant perturbations. One of these is the usual 0^ 
operator, whose coefficient renormalizes as. 



(3.9) 



Here, r always denotes the deviation from the critical point. The other relevant perturbation, 
whose discussion we have omitted thus far, is the chemical potential. 



SL = -nij-iij: 



laria 



(3.10) 



However, this perturbation is redundant, as it can be absorbed into a shift of hot spot 
location. Moreover, as already observed in section II, the location of the hot spots is not 
renormalized in the low-energy theory, which implies that there is no mixing between the 
two relevant operators. This is unlike the situation for the Ising-nematic transition in a 
metal studied in Ref. 16, where such mixing leads to a nontrivial shift of the Fermi surface 
as a function of deviation r from the critical point. 



Introducing the renormalized one-particle irreducible correlation functions of Uf fermion 
and Uf, boson fields 

we can write down the renormahzation group equations, 

(^|l + ^-a^ + ^"^ + ^-^^^ - ^^^^^ _ !^ _ !!|i^ r"-^ ({p}, a., u., 7. r., /.) = 

(3.12) 
Here, the /3-functions and anomalous dimensions are functions of Or and Uj. given by, 

dar 



(3a = fJ'- 



^ 1 7 

V^ = /^ o- log A 



dUr 

, Pu = Ai^— 

a,u,A Ofi 



a,u,A 

d 



1 d-fr 
7r OfJ' 



' '^^ = ^^TT ^°S Z^ 

a,u,A OH 



a,u,'y,A 

d. Zr 

, ^2 = /^7pl0g — 

a,u,A C/X /j(j) 



a,u,A 



(3.13) 
(3.14) 
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Using dimensional analysis, 

(3.15) 
Now, solving the RG equation (3.12), 

r-"/({cD},{p},a.,«.,f) = /-3"//2-"^Z^(.)-"^/^Z^(s)-"^-/^Z,(.)"^/^+"^/^-^ 

X r'""/(s-%(s){c:;},s-i{p},a,(s),M,(s),Z,(s)f) (3.16) 



with 



S—^ = /3a(ar(s),Ur(s)), a,.(l) = «,., S— ^ = /5„(ar(s),Wr(s)), Ur(l)=Ur 

as as 

f f" ds' ~ , \ / r ds' 

Z^{s) = exp ( / —i-rj^{ar{s'),Ur{s')) \ , Z^{s) = exp ( / — ^?7^(a,.(s'), ^^(s')) 

f f" ds' ~ \ f r ds' 

Z^{s) = exp ( / ^?7^(a^(s'), Mr(s')) J > ^r(s) = exp ( - / —^ri2{ar{s'),Ur{s')) 



(3.17) 



Now, let us construct the scaling forms of the correlation functions assuming that the 
couplings a^, u^ have a stable fixed point. Actually, as we will see below, this assumption 
is not supported by explicit calculations of low loop contributions to the /3-functions and 
anomalous dimensions. However, as already remarked, higher loop diagrams, which are 
naively suppressed by powers of 1/A^, actually scale as progressively higher powers of A^ 
and might modify the RG flow significantly. Thus, the fixed-point form of the correlation 
functions satisfies. 

Hence, typical frequencies and momenta are related by a; ~ |j)|", with the dynamical critical 
exponent z being given by, 

z = 2-ri^ (3.19) 

Moreover, the correlation length ^ away from the critical point scales as ^ ~ r~^ with 

1 



(3.20) 



2 + r/2 
Specializing to boson and fermion two-point functions, 

D-\uj,p)^r^^-^^^K{u:C.m ^^ W-'^^kiu/W.P) (3.21) 

G-\uj,p)^r^^/^-^'^^L{uji^pi) ^^ W/^-^'^L{uj/\p\\p) (3.22) 
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FIG. 5. The leading contribution to the ferniion self-energy. 



Here, the expressions on the right give the correlation functions at the critical point to which 
we confine our attention from here on. From Eq. (3.22) we may infer the fate of the Fermi 
surface at the critical point. We expect that as ^ — ?■ oo the Fermi-surface remains sharply 
defined. Close to the hot spots, the Fermi surfaces of fermions ipi and ip2 will evolve into 
straight lines with a fixed angle between them. At the hot spot, the fermion self-energy 
takes the form, 

G-i (w, p = 0) ~ w V2-r,^A (323) 

which is generally non Fermi-liquid like. On the other hand, away from the hot spot, if we 
define p± as the distance to the Fermi surface and py as the distance to the hot spot, for 
P± ^ P\\ and a; ^ pjj, we expect well-defined Landau quasi-particles, 

G(uj,p) (3.24) 

zu — vfP± 

with the Fermi velocity v and quasiparticle residue Z vanishing as we approach the hot spot 



as, 



n^-l C?/'^ N ^^Z^+'^V 



vf{p\\) -^ pl-\ Z(p||)~pJ'^"^'^ (3.25) 

The remainder of this section will provide a computation of the 4 renormalization con- 
stants Z^, Z^, Z^, Zy to leading order in 1/A^. At this order, the constants will depend only 
upon the dimensionless constant a, and do not involve u. We discuss the renormalization 
of u in Appendix B 2. Thus our considerations here will involve the RG fiow only of the 
single coupling a, the ratio of the velocities, and a discussion of its physical implications. 
For completeness, we will also compute the renormalization constant Zr, which determines 
the scaling of the correlation length away from the critical point. This constant will depend 
upon both a and u already at leading order in 1/A^. 

As we will see below, the 4 renormalization constants will be overdetermined from the 
structure of the 1/A^ corrections to the fermion self energy, the boson- fermion vertex, and 
the boson self energy. Computations of these quantities are provided in the appendix, and 
we use the results here to compute the Z's. 

The first correction to the self-energy of the fermion ipi^^ is given by Fig. 5, and computed 
in Appendix A 2. 
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2 / 2 A \ 

Si(w,p) = - ^ ,,1^1 isgn(a;)(v/7|a;| + {v2 ■ p)^ - \v2 ■ p\) + -h ■ P^og - — - (3.26) 

Note that unless otherwise stated, we will discuss the £ = 1 hot spot and drop the index i. 
We see that at the hot spot, p = 0, the self-energy has a non- Fermi liquid forni,^°'^° 

This result is consistent with our scaling form (3.23); to this order the anomalous dimension 
rj^ = 0. On the other hand, away from the hot spot, in the regime 7|c<;| ^ {v2 ■ p)^, the 
fermion propagator takes the Fermi-liquid form (3.24). To leading order, the Fermi surface 
is given by Vi-p = 0. The Fermi velocity and quasiparticle residue vanish with the distance 
p\l along the Fermi-surface to the hot spot as, 

VF = ^P\\, Z = — (27™)^V/M- + «1 P\\ (3.28) 

consistent with the scaling form (3.25) with mean-field exponents z = 2, rj^ = 0. 

The last term in Eq. (3.26) contributes to the renormalization of Vx, Vy, and so constrains 
the renormalization constants by 

Z^Z: = 1 - -^^^log(A//i) (3.29) 

Z^Zy = 1 + -^^^log(A//i) (3.30) 

irnJM 1 + a^ 



Next we consider the correction to the boson-fermion vertex, 

-{^2Ap)4.'iP)ri-q))iPi = C,r^^^^t(p, q){27rf6%p' - p - q) (3.31) 

This is given by Fig. 6 and computed in Appendix A3. We need only the UV divergent 
part, which is 

r,,,,t(j>,g) = l + ^tan-ilogA (3.32) 

Eq. (3.32) constrains the renormalizations by 



nnN a 



Z'fZ^ = 1 - ZTT? tan-' - log(A//i) (3.33) 



Finally, we consider the corrections to the boson two-point function, shown in Fig. 7, and 
computed in Appendix A 4. These yield 
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FIG. 6. The leading correction to the boson-fermion vertex. 




(a) 



(b) 




(c) 



(d) 



FIG. 7. The leading correction to the boson polarization. A sum over both directions of the fermion 
loop is implied. 
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u I loffA 



(3.34) 



Note that both the frequency and momentum dependent parts of the boson propagator 
receive renormalization corrections. As we discussed earlier, the corrections to the coefficient 
of |a;| should not be considered as renormalizations of an independent coupling 7, but should 
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rather track the renormahzations of the fermion velocities. Consequently, from Eqs. (3.8) 
and (3.34), we conclude that 

Z^{Z:Ziy^ = 1-^7 *^^"' - log(A//i) (3.35) 

From the momentum dependent part of (3.34) we immediately obtain the bosonic field 
strength renormalization, 

while the r dependent part of (3.34) yields the renormalization constant Zr, 



We note that while our results for the fermion self-energy (3.26) and the vertex (3.32) 
are in agreement with Ref. 14, the expression for the boson two-point function Eq. (3.34) 
differs from that of Ref. 14. More precisely, the frequency dependent part of our D~^ 
agrees with Ref. 14, while the momentum dependent part does not. As already noted, the 
renormalization of the frequency dependent part of D~^ is constrained by that of the fermion 
self-energy and the vertex. On the other hand, the renormalization of the momentum 
dependent part is completely independent. The authors of Ref. 14 found that both the 
frequency and the momentum parts are renormalized by the same factor, which would 
imply that the dynamical critical exponent 2; = 2 to this order. However, our calculations 
indicate that the two renormahzations are equal only at a = 1 and, as we will see below, 
the dynamical critical exponent z receives corrections already at the present order in 1/A^. 



We now have 5 equations for 4 renormalization constants: Eqs. (3.29), (3.30), (3.33), 
(3.35), and (3.36). It is easily verified that they are consistent with each other. This is a 
strong check on our renormalization procedure, and verifies the consistency of tying 7 to the 
velocities by Eq. (2.19). We can solve these equations to obtain 
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A. RG flows 



The renormalization constants in Eq. (3.38) determine the flow of the dimensionless 
coupling a with the /3-function 

12 o? 
/3(a,) = ^^ "'^ ^ (3.39) 

The /3 function for the velocity anisotropy a has an infrared stable fixed point a = and an 
infrared unstable fixed point a = oo. Physically, both fixed points correspond to a nested 
Fermi surface. For a = 0, the Fermi-velocities at the two hot spots are anti-parallel, while 
for a = oo they are parallel. The flows to the two fixed points are logarithmic. In particular, 
near the infrared stable fixed point a = 0, 

«.(s) = yfz-^ (3.40) 

1 + — ^log(l/s) 

Here we've assumed that the starting point of the flow a,. ^ 1. Note that the logarithmic 
flow to a — )■ in the infrared, with vanishing velocity ratio, is similar to that found recently 
in Ref. 19 in a different physical context. 

Let us now discuss the physics of the a = fixed point. The renormalization constants 
in (3. 36), (3. 37), (3.38) also determine the renormalization of the velocities, the anomalous 
dimensions of the bosons, fermions and of the 0^ operator. For the velocities, the ratio is 
already specified by a, and it is convenient to take 7 as the other independent combination 
of the velocities. We have therefore 

'77 = 77 ( «r ) ( 1 + ( "r ) tan~^ — 

■KnN \ar J \ \ar J ar 

' — ttr + — ^ + tt.. tan — 



irnN \ar \a~ ^ J a^ 



r 



7]^ = — I a,. ] ( 1 + ( a^ ] tan ^ — 

Note that as can be seen from Eqs. (3. 16), (3. 19) the flow of the dimensionful constant 7^ 
described by the exponent rj^ is equivalent to an anomalous dynamical critical exponent 
z. Since 77^ is non-zero, the dynamical behaviour of the theory deviates from the simple 
Hertz- Millis scaling with z = 2. 
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FIG. 8. Modification of the Fermi surfaces in Fig. 2 at the SDW quantum critical point. As in 
Figs. 2 and 3, the full lines are the Fermi surfaces, and the white, light shaded, and dark shaded 
regions denote momenta where 0, 1, and 2 of the bands are occupied. The equation of one of the 
Fermi surfaces is given in (3.43). 



As a flows slowly to 0, the critical exponents in Eq. (3.41) slowly vary: 

^^^iv^' ''^^-^h; ^^^;^^' "^'^-^h; "^-^° ^'-^'^ 

Observe that the corrections to the critical exponents diverge as a^ — ?■ 0. Thus, for suffi- 
ciently small momenta the 1/A^ expansion breaks down. From Eq. (3.42) we see that this 
will happen when a ~ 1/yN] from Eq. (3.40), we can estimate that this occurs at a mo- 
mentum scale k ~ exp(— A^^/^). This is parametrically smaller than the scale k ~ exp(— A^) 
at which the direct expansion in 1/A^ (without RG improvement) becomes invalid. 

Despite the breakdown of the RG at the longest scales, there is an intermediate asymptotic 
regime, 1/yN <C ttr ^ 1, where Eq. (3.42) remains valid, and we can integrate the RG 
equations and find interesting consequences for both the fermionic and bosonic spectra. 

For the fermions, the location of the ipi Fermi surface is given at tree-level by -Oi ■ p = 0, 
or Py = —v^px/vy = —px/a. Evaluating a at s = fi/px, we find the Fermi surface at 

12 
Py = - -:i^P^ log(/^/ \Px I ) (3.43) 

The resulting Fermi surface distorts from the shape shown in Fig. 1 to that in Fig. 8. We 
may also use RG to improve the one-loop result for the fermion self-energy (3.26). From 
Eq. (3.16), the fermion self-energy at the hot spot is, 

E(a;,p = 0) ~ -^exp (—^^\og'J^\ \uj\"\gn{uj), (3.44) 

Along the Fermi surface away from the hot spot, the quasiparticle residue and Fermi velocity 
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behave as, 

1/9 

The characteristic frequency of the bosonic spectrum is cj ~ (f /ir'-, evaluating 7,. at 
s = fi/\q\, we find that it scales with a 'super power-law' of the momentum 



-~«^^exp^j5^.og^|J. (3.46) 

From Eq. (3.16) we also obtain the static and dynamic scaling of the bosonic propagator, 

D-\u = 0,g1 ~ Igfexp (^^ log^ |) (3-47) 

Note that the unusual super-power law dependencies in Eqs. (3.44), (3. 45), (3. 46), (3. 47) 
are consequences of the scaling of a,. — )■ in the infrared and associated divergences of the 
anomalous dimensions. 

IV. COUNTING POWERS OF N 

As written in Eq. (2.17), our field theory offers a potentially simple way of organizing 
perturbation theory in powers of 1/A^: each boson propagator comes with a power of 1/N, 
each fermion loop yields a power of A^, and each u interaction yields a factor A^: we refer to 
this as the "naive" 1/A^ expansion, and it has been the basis of our computations so far. 

However, because we have to take 77 — )■ in the scaling limit, there is a danger that some 
of the higher order diagrams will have a singular dependence on 77. The fermion propagators 
in such diagrams need to include self-energy corrections for the diagrams to be finite in the 
?7 — !■ limit. The price we will pay for this regularization is that the diagram will acquire 
additional powers of A^, and the naive counting of powers of 1/A^ will break down. 

Recently, in the context of a theory of a Fermi surface interacting with a gauge field, 
S.-S. Lee^^ has given a procedure for identifying diagrams with a breakdown of naive 1/A^ 
counting, and shown that the expansion in powers of 1/A^ is actually an expansion in the 
genus of a surface defined by the graph. Using his methods we will show that many similar 
issues appear in our theory for the SDW transition of a Fermi surface, although subtle 
differences in RG properties imply that in the present case no genus expansion exists, and 
diagrams of increasingly higher order in A^ are generated as the number of loops is increased. 

In the absence of an external pairing vertex (see section V), the simplest diagrams ex- 
hibiting the above effect are the three-loop corrections to the boson-fermion vertex, see 
Fig. 9. In fact, the two diagrams are equal as they are related by particle-hole symmetry. 
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FIG. 9. Three loop corrections to the boson-fermion vertex that are enhanced in A^, scahng as 



The external fermions are taken to have hot spot index i = 1, while the fermions running in 
the loop can come from any hot spot i', although we will see that the singular contributions 
will originate from f = 1 and i' = 3. The diagram is given by, 



<^r,^.^t(p,g)r" = -r-r-r 






(2vr; 



6 



Gi{k)G2{k')D{k' - p)D{k - k')D{p + q-k) 
Substituting the four-point boson vertex /, Eq. (2.15), 



'^r^^,^t(p,g) 



0-02 V"! 



IN sr^ f dkrdkdk'dk' ,. . . ,,. .,, , . ., ., 



27rV^Vy 



(27r) 



{ir]{pr - kr) - vf ■ {p- k)){irj{qr + Pr - k'^) - t7f ■ (g + p- k')) 

X ^^ ^^D{k' - p)D{k - k')D{p + q - k) (4.1) 

{irjkT- — Vi ■ k){irik'^ — v^ ■ k') 

Observe that if £' = 2 or £' = 4 the four denominators in Eq. (4.1) involve four linearly 
independent combinations of internal momenta k, k'. As a result, the integral has a well 
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defined limit when 77 — )■ 0. On the other hand, when i' = 1 or i' = 3 (which we will 
also denote as i' = —1), v^ and v are parallel. Keeping only these two hot spots, let us 
integrate over the momentum components Vi ■ k, V2 ■ k'. We focus on the contribution from 
the fermionic poles, which, as we will see, is infrared singular. 

JN ^-^ f dkT-dk\\dk'^dk'n 

mr) - eji'ikr - Pr))){9ik',) - g(f (fc^ - p. - g.))) 

"" (ir/((l - f )A;, - pr) + m ■ p)(^r7((l - i')k'^ - Pr - Qr) + 1% ■ {p + q)) 
D{k' - p)D{k - k')D{p + q-k). 

Here k\\, k',, denote the components of k, k' along the Fermi surface of ipi and ip2 respectively, 
and the arguments of boson propagators are evaluated at vi ■ k = V2 ■ k' = 0. (Strictly 
speaking, only one pair of poles has Vi- k = V2- k' = 0, while the other pair has Vi-k = Vi-p 
and V2 ■ k' = V2 ■ {p + q). However, in situations of interest to us discussed below the above 
difference may be neglected in the bosonic propagators). 

Note that if we take the initial and final fermion momenta to lie on the Fermi surface, 
i.e. Vi ■ p = 0, V2 ■ {p + q) = 0, then ST diverges as 77"^. Since the dimension of 77 is w"^/^, 
this is synonymous to an infra-red divergence, 

S^,^M ^ V-'N-'u-\ (4.2) 

This behavior can be easily checked by, for instance, setting all the external momenta to 
zero {i.e. taking the external fermions to be at the hot spots). We also note that in the case 
when the external fermion momenta do not lie on the Fermi surface, the limit 77 — )■ can be 
taken in the contribution of hot spot pair i' = 1, but not i' = —1, as the latter contains a 
non-local UV divergence. Keeping 77 finite, we obtain, 

^^,^M ~ ^"'^"'Pl'- (4-3) 

where p± schematically denotes the distance of external fermion momenta to the Fermi 
surface. 

The infra-red divergences in Eqs. (4.2), (4.3) are a product of the bare fermion propagator 
having z = \ dynamics, whereas we expect that the full fermion propagator has the same 
dynamics as the spin-density wave excitations. We saw that this, indeed, holds at the one- 
loop level, where both the boson (3.2) and fermion (3.26) propagators are invariant under 
scaling with z = 2 (up to logarithmic corrections in the latter case). As in Ref. 15, the 
divergence can be cured by including the one-loop fermion self-energy within the fermion 
propagators, before taking the 7/ — )■ limit. This is the approach that will be adopted 
below. From Eq. (3.27), we know that the self-energy is ~ y/uJ/N. Therefore, mapping 
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rju — )■ ^/uj/N, we find from Eq. (4.2) tliat 

'^r.^^.^t ~ o(i) (4.4) 

Thus, tlie vertex correction is not suppressed relative to tlie bare value, and the naive 1/A^ 
expansion has broken down. In the appendix B 1, we compute the vertex correction in Fig. 9 
with dressed fermion propagators and find to logarithmic accuracy, 

^r,^.^I~^(«)i°g^ (4-5) 

where X is a finite negative function of a. Note that the strong infra-red divergence of 
Eq. (4.2) is now replaced by a mild logarithmic divergence that one may hope to treat with 
renormalization group. However, the price one has to pay for curing the strong infra-red 
divergence is the enhancement of the diagram with A^, as anticipated in Eq. (4.4). This 
enhancement occurs for any external fermion momenta (not only for momenta on the Fermi 
surface). Finally, the presence of a logarithm implies that not only is the diagram itself 
unsuppressed relative to its bare value, but also that the anomalous dimensions are not 
expected to be suppressed with A^. 



Having seen an explicit example of violation of naive large- A^ counting, we would like to 
investigate the general scaling of diagrams with A^ in our theory, when a one-loop dressed 
fermion propagator is used. Our procedure closely follows that of Ref. 15. A general diagram 
can be schematically written as, 

v=N^f ffi d\du, n - — ij^—^T n ^(^'^) (4-6) 



i = l j: 



Here, // and lb are numbers of fermion and boson propagators respectively, Lf is the number 
of fermion loops and L is the number of total loops. The momenta Ij and qk are linear 
combinations of pi entering the fermion and boson propagators. The "naive" scaling of the 
diagram with A^ is given by "D ~ A^*^", 

Qo = Lf - h (4.7) 

It is clear that the enhancement of diagrams with A^ comes from the dangerous factor of 1/A^ 
in the fermion self-energy. However, in order to access this factor the fermion momentum 
must be on the Fermi surface. Given a diagram, let us call the phase-space for all internal 
fermion momenta to lie on the Fermi surface, the "singular manifold." Having identified this 
manifold, one can divide the momentum integration variables into components parallel p\\ 
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and perpendicular p± to the manifold, 

L n 2L—n 

Yld^Pi = l[dp\\al[dp^b (4.8) 

j=l a=l 6=1 

where n is the dimension of the manifold. Linear combinations of p_l's enter the fermion 
energy v ■ Ij and hence scale as 1/A^, making the fermion propagators scale as A^. On the 
other hand, the components p\\ only enter the bosonic propagators and the one-loop fermion 
self-energy T^uoop and scale as N^. Hence, the diagram acquires an enhancement, V ~ A^*^, 
Q = Qo + AQ, 

AQ = [If -2L + n] (4.9) 

where [x] = x if x > and [x] = if x < 0. 

Thus, to find the degree of a diagram in A^, one has to find the singular manifold and 
compute its dimension n. This can be done diagramatically by introducing a double-line 
representation, originally used in the study of electron-phonon interactions.^'' Below, we will 
consider diagrams involving opposite hot spot pairs £ = 1 and £ = —1 only. Subsitution of 
fermions from hot spots £ = 2 and £ = —2 into these diagrams is expected to reduce the 
dimension of the singular manifold. Moreover, we for simplicity consider diagrams without 
the quartic bosonic vertex u. Finally, we take all the exernal fermion momenta to be on the 
Fermi surface. 

Now, we are ready to introduce the double-line representation. We would like to find 
under what conditions do all the fermions in a diagram go to the Fermi surface. Observe, 
that any momentum can be uniquely decomposed into components along the Fermi surface 
of fermion 1 and fermion 2. Thus, we fatten bosonic propagators into double lines, one 
carrying momentum along the Fermi surface of fermion 1, and the other along the Fermi 
surface of fermion 2. If a fermion is to absorb this bosonic momentum and stay on the Fermi 
surface, its incoming and outgoing momenta are fixed in terms of the components of the 
double line. Hence, the boson-fermion vertices can be redrawn as shown in Fig. 10. Note 
that if a certain momentum is along the Fermi surface of fermion 1 from hot spot £ = 1, 
it is also along the Fermi surface of fermion 1 from hot spot £ = —1. Thus, the fermion 
lines in our diagrams can come from either of these hot spots. Also, the direction of lines 
in the double-line representation is not fixed, and need not coincide with that in the single 
line representation. If the two are opposite, then it is understood that the physical fermion 
momentum p is the negative of the momentum carried by the fermion in the double-line 
representation, see Fig. 11. Because we are neglecting the Fermi surface curvature in the 
low-energy theory, a particle with momentum p is on the Fermi surface if and and only 
if a particle with momentum —p is on the Fermi surface, and the above representation is 
consistent. (We remind the reader that here all the fermion momenta are defined relative to 
hot spot locations). 
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1 ^ 1 

FIG. 10. Double line representation for the boson-fermion vertex. 



P= -P2 




FIG. 11. Double line representation for the boson-fermion vertex. The direction of momentum 
and particle flow need not coincide. 



Thus, the double line representation completely specifies the singular manifold. In partic- 
ular, the dimension of the manifold n is just given by the number of loops in this representa- 
tion. As an example, consider the double line represenation of the diagrams in Fig. 9 shown 
in Fig. 12. We see that Fig. 12 contains two closed loops, which implies that the singular 
manifold is two-dimensional. From Eq. (4.9), the enhancement of the diagram is AQ = 2, 
which combined with the naive degree of the diagram, Eq. (4.7), Qq = —2, gives Q = 0, 
consistent with the explicit calculation in Eq. (4.5). In Fig. 13 we also give an example of a 
vertex correction which is not enhanced in A^. Here, the double line representation contains 
no loops so the dimension of the singular manifold is zero, AQ = and the degree of the 
diagram is given by the naive A^ counting, Q = —2. 

It is easy to see that the violations of naive large- A^ counting are not confined to vertex 
corrections alone. In Fig. 14 we show a fermion self-energy diagram that acquires an en- 
hancement. Indeed, the naive degree of the graph is Qq = —3. However, since the double 
line representation contains three loops, the graph receives an enhancement AQ = 2, so that 
the total degree of the graph is Q = —1. Hence, the graph is of the same order 0{1/N) as 
the one-loop fermion self-energy. Similarly, in Fig. 15 we show an enhanced diagram for the 
boson self-energy. In this case, Qq = —1, AQ = 2, Q = 1. Hence, the diagram is of 0{N), 
again the same as the tree level contribution. 

A remarkable feature of the large- A^ counting in Eqs. (4.7), (4.9), pointed out in Ref. 15, 
is that the degree of a diagram is related to its topology. Let us first apply the topological 
classification to vacuum energy diagrams, i.e. graphs with no external lines. We can convert 
these diagrams into two-dimensional surfaces in the following way. First, let us introduce 
fermion loops back into the double line represenation (they will appear dotted in our dia- 
grams, see Fig. 16). Then attach a face to each solid loop of the double-line representation 

25 




FIG. 12. Double line representation applied to the diagrams in Fig. 9. The enhancement of the 
diagram in N is related to the number of loops n in the double line-representation via Eq. 4.9. 




FIG. 13. A three loop vertex correction with no enhancement in A^. 



and a face to each dotted loop {i.e. ferinion loop). As a result, each boson propagator is 
shared by two faces with solid boundaries, while each fermion propagator is shared by a face 
with a solid boundary and a face with a dotted boundary. Therefore, if we glue the faces 
along propagators we obtain a closed surface. Now consider the Euler characteristic of this 
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FIG. 14. A diagram for the ferniion self-energy that is of 0{\/N) as a result of enhancement. 




FIG. 15. A diagram for the boson self-energy that is of 0{N) as a result of enhancement. 



surface, 



X=F-E+V 



(4.10) 



where F is the number of faces, E is the number of edges and V is the number of vertices 
of the surface. We have, F = Lj + n, E = I^ + If and V is just the number of vertices in 
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FIG. 16. Converting vacuum energy diagrams into surfaces: a face is attached to each sohd and 
grey loop in the double-hne representation (on the right). In the present case, the resulting surface 
is a sphere. 



the original Feynman graph. Now, using V = 2Ib, 2V = 21 f we obtain, 

V 

X = Lf + n-- 



(4.11) 



However, using L = Ii, + If — V + l,we see that the degree of a diagram in N given by Eqs. 
(4.7), (4.9), is, 

(4.12) 



V 
Q = Lf-- + n-2 



where we've assumed that the argument of [] in Eq. (4.9) is positive. Thus, we arrive at the 
relation, 

Q = X-2 (4.13) 

This result means that at each order in 1/A^ one has to sum an infinite set of diagrams with 
a given Euler characteristic. In particular, at A^ = oo the theory is dominated by diagrams 
with X = 2, i.e. those whose double-line representation can be drawn on a sphere. Such 
graphs are often referred to as planar diagrams. 

It is straightforward to extend the classification above to diagrams with external legs. For 
instance, fermion self-energy diagrams can be obtained by cutting one fermion propagator 
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in a vacuum graph. This results in lb ^ lb, Lf ^ Lf ~1, so Qo ^ Qq — I, and // — )■ J/ — 1, 
L— T-L — 1, n ^ n — 1, as cutting a fermion propagator destroys a sohd loop in the double 
line representation. Hence, AQ — )■ AQ and Q — > Q — 1, «.e. 

g = X-3 (4.14) 

with X the Euler characteristic of the initial vacuum diagram. In particular, planar vacuum 
graphs give rise to fermion self-energy diagrams of 0{1/N). 

Similarly, to obtain a boson self-energy diagram, we cut a boson propagator in a vacuum 
bubble. This gives lb ^ h — 1, Lf — )> Lf, so Qo ~^ Qo + 1; and If ^>- If, L -^ L — 1, 
n — )■ n — 2, as we now destroy two solid loops in the double line representation. Hence, 
AQ ^ AQ and Q ^ Q + 1, i.e. 

Q = X-1 (4.15) 

Hence, planar graphs give rise to boson self-energy diagrams of 0{N). 

Likewise, to obtain vertex correction diagrams, we remove a vertex in a vacuum bubble. 
As a result, /{, — )■ /{, — 1, Lj — )• Lj — 1, so Qq — ?■ Qo, and Jj — )■ Jj — 2, L — )■ L — 2, n — )■ ri — 2, 
as we again destroy two solid loops in the double hne representation. Hence, AQ — )■ AQ 
and Q ^ Q, i.e. 

Q = X-2 (4.16) 

and all planar graphs give rise to vertex diagrams of (9(1). 

At this point, we would like to make a remark about conditions on external momenta in 
diagrams needed for the enhancements to occur. Up to now we have been assuming that all 
the external fermion momenta in a diagram are on the Fermi surface. If all the diagrams in 
our theory were UV finite then this condition would, indeed, be required. However, as we 
have seen, some of the diagrams actually contain logarithmic divergences, i.e. they receive 
contributions from momenta, which are much larger than the external momenta. For the 
purpose of computing the UV divergent contribution to these diagrams and estimating its 
scaling with A^, we can set the external momenta to zero (which certainly puts the external 
fermions on the Fermi surface). This explains why the vertex correction in Figs. 9,12 receives 
an enhancement for any external fermion momentum, as can be explicitly seen in Eq. (4.5). 

So far, we have left out one type of diagram which is important from the point of view 
of RG properties of the theory, namely diagrams for the boson four-point function. Such 
diagrams can be obtained by cutting two boson propagators in a vacuum bubble. This 
results in Ij, ^ h — 2, Lf -^ Lf, so Qo — ^ Qo + 2. Now let us discuss the change in the 
enhancement AQ. We see that // -^If, L — )■ L — 2. The change in the dimension of the 
singular manifold Sn depends on how many loops in the double line representation the two 
propagators that we cut share. If both the components 1 and 2 of the two propagators 
are part of the same two solid loops, see Fig. 17c, then the change in the dimension of the 
singular manifold 6n = —2. If these two propagators share only one solid loop, see Fig. 17b, 
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(b) 





(c) 




FIG. 17. Producing a boson four-point function from a vacuum bubble by cutting two boson 
propagators. If the initial diagram is planar and only two solid lines are cut in the double-line 
representation then the resulting diagram is disconnected, as in (a). Diagrams of highest degree 
are obtained by starting with a planar diagram and cutting three solid line loops, as in (b), or 
starting with a diagram with x = 1 and cutting two solid line loops, as in (c). 



then 5n = —3. Finally, if the two propagators don't share any solid loops, then Sn = —4. 
Thus, we obtain, AQ — )• AQ + 4 + 6n and Q ^ Q + 6 + 6n, i.e. 



Q = X + 4 + 6n 



(4.17) 



It appears that the highest possible degree of the four-point vertex corresponds to starting 
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FIG. 18. A diagram for the boson four-point function that diverges logarithmically and scales as 



with a planar graph and cutting two bosonic propagators, which are part of the same double- 
line loop, to obtain, Q = 4. However, it is easy to see that this always produces a diagram, 
which is disconnected, see Fig. 17a. To obtain a connected diagram for the four-point 
function starting from a planar graph, we must cut at least three solid loops, such that 
the highest possible degree of a four-point function is Q = 3. The fact that the four-point 
vertex scales as N^ could be anticipated from the simple one- loop result in Eq. (2.15). 
Indeed, for special kinematic conditions, vi ■ (g2 + ^3) = 0, t?2 ■ (^1 + ^) =0, Eq. (2.15) 
diverges as Nirfu)'^ , which after including the one- loop fermion self-energy is expected to 
become of order A^^. Such kinematic conditions are automatically assumed in our double 
line representation that led to the large- A^ counting in Eq. (4.17). However, as was already 
noted, diagrams that have ultraviolet divergences are expected to receive the enhancement in 
Eq. (4.9) independent of external momenta. The simplest diagram for the boson four-point 
vertex that is expected to scale as N^ and exhibits such a divergence is shown in Fig. 18. 
In the appendix, we explicitly evaluate this diagram obtaining to logarithmic accuracy. 






(4.18) 



with Y a finite function of a. 

The fact that there are diagrams for the four-point boson function that scale as A^^ for 
arbitrary external momenta has drastic consequences for the theory. Indeed, a diagram with 
just quartic internal vertices (which can themselves have a non-trivial internal structure). 
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will scale as A^*^ , with Q = V4, + ^, where V4 is the number of quartic vertices and Eb is the 
number of external bosons. Thus, the degree of the diagram in A^ grows with the number of 
quartic vertices. This means that perturbation theory based on the one-loop dressed fermion 
propagator is not a good starting point for taking the large- A^ limit, and no genus expansion 
similar to that of Ref. 15 exists in the present case. Note that this effect was not captured 
in our initial large- iV counting, as we have ignored the possible presence of UV divergent 
sub diagrams. 



V. PAIRING VERTEX 

In this section we will study the renormalization properties of the BCS order parameter 
to one loop order. We consider pairing in the spin singlet, parity even, momentum zero 
channel. There are four order parameters that one can form out of our four pairs of hot 
spots. 



K. = e..'(^f:-vf=/ + /xC-v^=/) + ^6..'(^fr'^f:'' + /^^^r'C'') (5.1) 



111/ 



Here the minus sign in the hot spot labels £=—1=3 and £=—2 = 4 denotes the opposite 
hot spot pair. The geometry of the pairing operators for £ = 1 is illustrated in Fig. 19. The 




FIG. 19. Pairing of the electrons at the i = ±1 hotspots of Fig. 1. Electrons at opposite ends of 
the arrows form spin-singlet pairs. The fj, = +1 {p, = —1) pairing amplitude in Eq. (5.1) has the 
same (opposite) sign on the two arrows. Only the /i = — 1 spin singlet pairing is enhanced near 
the SDW critical point. 



coefficients /i = ±l, z/ = ±l determine the transformation properties of V under the lattice 
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TABLE I. Symmetry properties of the pairing vertex. 
rotation symmetry -R^/2 and the reflection symmetry /(-i,i) about the (—1, 1) axis: 



R-k/2 ■ V^u -^ VV, 



fii/ 



-^(-1,1) : V^^ -)■ /iV^j, 



(5.2) 
(5.3) 



These properties are summarized in Table I. Since the theory (2.17) conserves the number 
of fermions at each hot spot pair i, the parts of the order parameter involving £ = ±1 and 
i = ±2 renormalize independently. Hence, the scaling dimension of the pairing vertex in 
the low-energy theory is independent of u and is sensitive only to /i, i.e the operators with 
s and dxy, and g and d^2_y2 symmetries are degenerate. 

The renormalization properties of the operator V can be determined from its insertion 
into the correlation function, 



e^^'J 



^D 



\i=^\. 



'M='^i 



rv'^tv,t(fci,A;_i) = d xid x_i{V{0)ipl^, (x_i)^|^ (a;i))ip/e 



J{kixi+k-ix-i) 



(5.4) 



At tree level, Fy^t^t = 1- Let us now consider the one- loop renormalization of V, shown in 
Fig. 20 a). This diagram is given by 



'5Fy^t^t(A;i, A;_i) 



-3yU 



dH 



■D{l)Gl{h - l)G^\k^i + I). 



(5.5) 



Details of the evaluation of (5.5) appear in Appendix B3. Direct computation with bare 
fermion propagators gives rise to strong infra-red divergences, which are cured by using the 
one-loop dressed propagators. With this approach, we obtain to logarithmic accuracy 



6r 



yi/)t^t 



/ia 



7r(a^ + Ij \1^ J 



(5.6) 



Note that the one loop renormalization of the pairing vertex (5.6) is of order unity, and is 
not suppressed in 1/A^. Thus the naive counting in powers of 1/A^ is violated, as was already 
noted in Ref. 12. Moreover, the one-loop contribution gives a suppression of the vertex for 
/i = 1 (s and dxy channels) and an enhancement for /i = — 1 {dx2_y2, g channels) as expected. 
Finally, we find that the one-loop result has a non-local log divergence. The origin of this 
non-local divergence is BCS pairing of the Fermi surface away from the hot spots. Indeed, 
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FIG. 20. The leading corrections to (a) BCS pairing vertex, (b) density-wave vertex. 



as noted in Appendix B3, the divergence comes from the regime where ^\It\ <^ ^y, with l\\ 
the component of / along the Fermi surface of ^2- This is precisely the regime in which one 
has good Landau-quasiparticles, suggesting that it may be possible to obtain Eq. (5.6) in a 
Fermi liquid computation. 

We now show this is indeed the case, and obtain (5.6) in a physically transparent form. 
Let us approximate the propagators in Eq. (5.5) by the Fermi-liquid form Eq. (3.24), 



oFy^t^t 



3^ ^di^ r du /"^^ ^W ^(/ji) 

iVy 2^y^|;^|<;2 2^7 ^l\lr\ + /J iilr - w) - vf{Ii\)1± i{lr + uj) + vf{1\\)Ii. 



(5.7) 
with the Fermi-liquid parameters given by Eq. (3.28). Note that due to the restriction 
7|/^| <^ /n the bosonic propagator is static. Changing variables to e = Vf{1\\)1±-, 

_3/i [dk z\k) r dW f de 1 1 . . 

^'^^'^^ N J 27rt;^(/||)/f^^|<,|27ry 27r t{l, - u) - e z{lr + uj) + e ^ ' ' 

The integral over l^-, e has the form familiar from Fermi-liquid theory and gives the usual 
BCS logarithm, 

dlr f de 1 1 1 ApL ir- nN 

log (5.9) 



l-K } 2n i{lr — u)) — e i{lr + Lo) + e 2n cu 
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where KpL is the frequency /energy cut-off, which in the present case is KpL = ^f\/l- Of 
course, for the above form to hold, we need u ^ A^^,. Thus, 

2'K^N J^^ "vf{Iii)Iii luj 7r(a2 + i) ^o; 

which agrees with the result in Eq. (B26) obtained from a more complete computation. Note 
that the prefactor of 1/A^ arising from the boson propagator has disappeared from the final 
result. A similar log-squared term has been noted for the pairing vertex in a theory of a 
Fermi surface coupled to a gauge field in three dimensions'^ '^^ and in a theory of a Fermi 
surface interacting via a Chern-Simons gauge field and a 1/r potential in two dimensions.'^ 

The appearance of the log-squared term above indicates a breakdown of the present RG 
in analyzing the renormalization of the pairing vertex. It is clearly a consequence of two 
different physical effects. One is the familiar BCS logarithm of Fermi liquid theory, which 
appears here from the Fermi surface away from the hot spots. The second logarithm is a 
critical singularity associated with SDW fluctuations at the hot spot. Our RG approach, 
defined in terms of a cutoff A which measures distance from the hot spot, is unable to 
regulate the first logarithm: the Fermi surface is present at momenta all the way upto A. 

An alternative RG is necessary to analyze the consequences of the log-squared term. One 
possible approach is that of Son'^, who worked with an RG defined in terms of momentum 
shells a fixed distance from the Fermi surface of fermions coupled to a gauge field. We leave 
such investigations for future work. 



VI. DENSITY VERTICES 

In this section we focus attention on one of the interesting consequences of the pseudospin 
symmetries of the critical theory of the SDW transition, specified by Eq. (2.6). Note that 
the pseudospin rotations can be performed independently on different pairs of hotspots. 

Under the operation in Eq. (2.6), the pairing operator (5.1) in the particle-particle channel 
becomes exactly degenerate with certain operators in the particle-hole channel which connect 
opposite patches of the Fermi surface. Indeed, consider spin-singlet operators that can be 
built out of fermions coming from hot spots i and —i. Using the spinor representation (2.3), 
we may write these as, 

V^^ = M,,e.„,%^^¥^^„, (6.1) 

The indices a, /3 of Vap carry spin 1/2 under the independent SU^^{2) and SU^{2) particle- 
hole symmetries. Hence, we have a set of four degenerate operators. Choosing a = 1, 

/5 = 1, 

V(, = M,,e..,C.V,v (6.2) 
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The mixing matrix Mij is fixed by lattice symmetries to give operators, 



T/£,Q=(0,0) . 

V 

T/£,Q=(T,7r) 

V 



Ccrcr' 






(6.3) 
(6.4) 



which correspond to superconducting order parameters with momenta (0, 0) and (vr, tt) re- 
spectively. The index // = ±1 determines the parity of the operator under a reflection 
about a lattice diagonal. The operator (6.3) was considered above. We will not discuss the 
other operator (6.4) below; due to kinematics, its renormalization at one-loop order contains 
neither the large- A^ enhancement, nor the unusual powers of logarithm squared. 

Now, let us discuss the particle-hole partners of (6.3). Setting a = 2, /3 = 2 in (6.1) 
simply gives rise to the Hermitian conjugate of (6.3). On the other hand a = 2, /3 = 1 gives 
the operators, 

(6.5) 



oi = r.^'i^l + ^^r.^Hl 



The other choice a = 1, /3 = 2 generates the Hermitian conjugates of (6.5). Following Fig. 19, 
the O^ operators are illustrated in Fig. 21. To determine the wavevectors of these operators. 




FIG. 21. Spin singlet density operators (~ ilvip) of the electrons at the £ = ±1 hotspots of Fig. 1 
(see also Fig. 19), shown with an arrow pointing from the Brillouin zone location of V' to that of 
ip. The dashed arrows are the density operators in the first Brillouin zone. The full arrows are in 
an extended zone scheme which shows that these operators have net momentum Qi = 2Ky{—l, 1), 
where [Kx,Ky) is the location of the ^ = 1, i = 1 hot spot. The density operator with opposite 
signs (/i = —1) on the two arrows is enhanced near the SDW critical point. Similarly the i. = ±2 
hot spots contribute density operators at Q2 = 2Ky{l, 1). 
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let the £ = 1, i = 1 hot spot be at Ki = {K^, Ky). (Note that here we are using the principal 
axes of the square lattice for the momentum co-ordinates, not the diagonal axes indicated 
in Fig. 1.) Then, from Fig. 1 we note that the £ = 1, i = 2 hot spot is at {—Ky^—K^), 
and so the value of the SDW wavevector Q = (tt, vr) implies that K^ + Ky = tt. Also from 
Fig. 1, the £ = —1, i = 1 hot spot is at {—Kx, ~Ky), and so we conclude that the ordering 
wavevector of the first term in Oj^ is (2Kj:,2Ky). Similarly, the ordering wavevector of the 
second term in O^ is seen to be {—2Ky, —2Kx). Using K^ + Ky = n, we observe that these 
two ordering wavevectors are actually equal, and take the common value Qi = 2Ky{—l, 1), 
which is therefore the momentum of the Oj^ order parameters, as shown in Fig. 19. Similarly, 
the momentum of the O^ order parameters is seen to be Q2 = 2Ky{—l, —1). Thus the Oj^ 
represent density modulations along the diagonals of the square lattice. 

For a clearer physical interpretation of the O^ orders, it is useful to express them in terms 
of the lattice fermions cg^, where the momentum k ranges over the full square lattice Brillouin 
zone. Then by looking at the transformations of Eq. (6.5) under all square lattice space group 
operations, and under time-reversal, we find that the O^ are orders are characterized by 

c^- - , Cr^n /9 ) = Oifo(k), (6.6) 

where /o(A;) is any periodic function on the Brillouin zone that is invariant under the point 
group operations which leave the wavevector Qi invariant i.e. under the little group of 
Qi. Also time-reversal and inversion symmetries imply /o(/c) is real and even. The little 
group consists only of reflections along the diagonals, and so a simple choice is /o(fc) = 
1 + Ci {co^kx + cos ky) + . . ., where Ci is a constant. By taking a Fourier transform of 
Eq. (6.6), it is clear that Of corresponds to an ordinary charge density wave (CDW) on the 
sites of the square lattice: 

£=1,2 

As we saw in Section V, SDW fluctuations suppress pairing with /i = +1, and so its particle- 
hole partner, the CDW order parameter O^ will also be suppressed. We will therefore not 
consider it further. 

By the same reasoning, the order parameter Oi should be enhanced by the SDW fluctu- 
ations, and so it is of far greater interest. Following the steps leading to Eq. (6.6), we now 
find 

(4-4/2,/fc+4/2,a) = O- /o(^) (Cosfcx - COS ky) , (6.8) 

where fo{k) has the same structure as foik). Time- reversal symmetry played an important 
role in constraining the rhs: it is easily verified that Eq. (6.8) is invariant under time-reversal 
for general complex 0^_. The order in Eq. (6.8) is odd under reflections along the diagonals, 
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FIG. 22. (Color online) Plot of the bond density modulations in Eq. (6.10). The lines are the 
links of the underlying square lattice. Each link contains a colored square representing the value 
of (cl^cgfj), where r and s are the sites at the ends of the link. We chose the ordering wavevector 

Qi = (27r/16)(l, —1). Notice the local Ising-nematic ordering, and the longer wavelength sinusoidal 
envelope along the diagonal. 



and so it is a px±y-densitj wave, in the nomenclature of Ref. 25. Despite the d- wave-like 
factor on the rhs of Eq. (6.8), this order is not the popular (i-density wave^^; the latter is 
odd under time-reversal, and in the present notation takes the form 



'^k-Q/2,a'^k+Q/2,c 



~ i (sin kx — sin ky) 



(6.9) 



with Q = (tt, tt). The order in Eq. (6.9) is not enhanced near the SDW critical point, while 
that in Eq. (6.8) is. By taking the Fourier transform of Eq. (6.8), it is easy to see that 
Oi does not lead to any modulations in the site charge density ( ct^c^o-), and so it is not a 

CDW. The non-zero modulations occur in the off-site correlations ( ci^cga ) with f ^ s. For 
r* and s nearest-neighbors, we have 



cA = Y^ fote^'^'-'-'^+'^/' + c.c.) [6. 



£=1,2 



■') — s,x ' ^s—r,x "r—s,y "s—r,y\ 



y\ ■> 



(6.10) 



where x and y are unit vectors corresponding to the sides of the square lattice unit cell. The 
modulations in the nearest neighbor bond variables (cf^^cp+x,^) and \c^„Cf^+y,a) are plotted 
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FIG. 23. (Color online) As in Fig. 22, but for orderings along both Qi = (27r/16)(l, — 1) and 
Q2 = (2^/16)(l,l). 



in Figs. 22 and 23. These observables measure spin-singlet correlations across a link: if there 
are 2 electrons on the 2 sites of a link, this observable takes different values depending upon 
whether the electrons are in a spin singlet or a spin triplet state. Thus Oi has the character 
of a valence bond solid (VBS) order parameter. The first factor on the rhs of Eq. (6.10) shows 
that the VBS order has modulations at the wavevectors Q^ along the square lattice diagonals. 
However, from our discussion above, note that \Qi\ = 2^/2Ky, where the magnitude of Ky 
is quite small for the Fermi surface in Fig. 1: the £ = 1, i = 1 hot spot is at {Kx,Ky). 
Thus the first factor in Eq. (6.10) contributes a relatively long- wavelength modulation, as is 
evident from Figs. 22 and 23. This long-wavelength modulation serves as an envelope to the 
oscillations given by the second factor in Eq. (6.10). The latter indicates indicates that the 
bond order has opposite signs on the x and y directed bonds: this short distance behavior 
corresponds locally to an Ising-nematic order, which is also evident in Figs. 22 and 23. The 
ordering in Eq. (6.10) becomes global Ising-nematic order in the limit Q^ — )■ 0. Non-linear 
terms in the effective action for the bond order will lock in commensurate values of Qi, and 
so it is possible that strong-coupling effects will prefer Qe = 0. 

As already remarked, the particle-hole symmetry of our theory guarantees a degeneracy 
between the d-wave superconducting vertex and the density-wave vertex. However, this 
degeneracy is lifted once effects which break the particle-hole symmetry are introduced. One 
such effect is the curvature of the Fermi surface at the hot spots. Nominally, the curvature 
is irrelevant under the scaling towards hot spots (2.16). However, we recall that the double- 
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log structure in Eq. (5.6) originates from an interplay between scaling in a Fermi-liquid 
and quantum critical scaling. Moreover, we know that the scaling of the superconducting 
vertex and the density- wave vertex in a Fermi liquid are very different: at one loop the 
corrections to former are logarithmic, while corrections to latter are suppressed by uj^l'^ . 
Thus, one might expect that the Fermi surface curvature will play an important role in the 
renormalization of the density-wave vertex, reducing its enhancement compared to the BCS 
vertex and establishing superconductivity as the dominant instability of the SDW critical 
point. We check this by an explicit calculation below. 

We introduce the Fermi-surface curvature into the theory via a perturbation, 

^^ = iEi(v-4^)^'i' (6-11) 

where rin ^ = z x -0^ is the unit tangent to the Fermi surface of ^/'f . 

Let us define the insertion of the density-wave order parameter O^ into the fermion 
correlation function. 



Vo^^^{k^.k^^)b,,, = J rf^xirf^a;_i(O^(0)^r;(x_iX,(xi))ip,e^^'^^^i-^-i^-i^ (6.12) 

At tree level rQ^^t(fci, fc-i) = 1. The one loop correction to the vertex is given by the 
diagram in Fig. 20b). We perform the calculations with propagators dressed by the one- 
loop fermion self-energy and by the curvature (6.11). Details are presented in Appendix B4. 
To leading logarithmic accuracy we obtain, 

^^ STT[a^ + 1) 70; 

which is a factor of 3 smaller than the corresponding expression for the superconducting 
vertex (5.6). 



Finally, we note the resemblance between our results and those obtained by Halboth and 
Metzner,^^ and Honerkamp et. al,^^ using a functional renormalization group treatment of 
the Hubbard model. They find dominant instabilities to SDW order and d-wave pairing, 
along with a sub-dominant enhancement of Ising-nematic order. They assumed their Ising- 
nematic order was at Q^ = 0, but their results could be limited by the finite resolution of 
Fermi surface points, and their specific Fermi surface configurations. It would be interesting 
if higher resolution studies of more generic Fermi surfaces lead to ordering compatible with 
Eq. (6.8). 
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VII. CONCLUSIONS 

Quantum phase transitions involving symmetry breaking in the presence of a Fermi sur- 
face can be associated with the appearance of a condensate of particle-hole pairs of the 
Fermi surface quasiparticles. Such transitions can be divided into two broad classes: those 
in which the particle-hole condensate carries net momentum Q 7^ 0, and those in which the 
particle-hole condensate is at Q = 0. Both classes were considered by Hertz in his 1976 
paper^, using a self-consistent RPA approach, formulated in terms of a RG analysis of an 
effective action for the condensate fluctuations. He argued that for both cases, and for all 
spatial dimensions d>2, the condensate fluctuations were effectively Gaussian, and hence 
the leading critical behavior could be exactly calculated. 

We have re-examined both classes of Fermi surface transitions in this and a previous 
paper^^. While Hertz's conclusions are expected to be largely correct in (i = 3, they break 
down^^ in both classes for the physically important case oi d = 2. Our previous paper^^ 
proposed and analyzed a critical theory in d = 2 for a paradigm of the Q = case: the 
onset of Ising-nematic order. This theory involved both the bosonic order parameter and 
the fermionic quasiparticles as fundamental degrees of freedom, which interact strongly at 
the quantum critical point. The present paper has considered a typical case in d = 2 with 
Q ^ 0, the onset of spin density wave (SDW) order, using a field theory for the bosonic 
order parameter and the fermions proposed by Abanov and Chubukov^°. 

Our analysis for Q ^ begins by focusing on the vicinity of the "hot spots" on the Fermi 
surface shown in Fig. 1. Zooming in on a single pair of hot spots, and shifting one of the hot 
spots by a momentum Q, we obtain the situation shown in Fig. 2, where we can approximate 
the two Fermi surfaces near the hot spots by two non-coUinear straight lines. The two Fermi 
surfaces are coupled at the hot spot by the SDW order parameter 0, and the low energy 
physics is then described by the field theory in Eq. (2.1). In the phase with SDW order 
with {(p) ^ 0, the Fermi surfaces reconnect into the configuration shown in Fig. 3, leading 
to electron and hole pockets appearing from the original large Fermi surface in Fig. 1. 

Our RG analysis of Eq. (2.1) was performed using the 1/A^ expansion, where the fermions 
are endowed with an additional flavor index which runs over A^ values. Initially, it seems 
that the counting of powers of 1/A^ is simple: each boson propagator comes with a factor 
of 1/A^, and each fermion loop yields a factor A^. Using this "naive" counting, all RG flow 
equations were computed to order 1/A^ in Section HI. We found a consistent renormalization 
of the couplings in the local field theory in Eq. (2.1); the damping parameter 7 appearing 
in the boson propagator was tied to the local couplings via Eq. (2.9), and this relation 
was maintained under the RG. The flow of the spin-damping rate under RG implies that 
the dynamical critical exponent z renormalizes away from its RPA value z = 2. This is 
in stark contrast to Hertz theory^ and previous studies of the present theory. ^^ One of the 
main consequences of the RG flow in Section HI was a logarithmic divergence in the ratio 
of Fermi velocity components with length scale: this implied that the Fermi surfaces at the 

41 



quantum critical point took the shape in Fig 8. The effective dynamical nesting of the Fermi 
surfaces at low energies gives rise to a divergence of anomalous dimensions, which may lead 
to a first order phase transition. 

Section IV looked at higher loop effects which showed that the naive counting of powers 
of 1/A^ was not correct. The enhancements in powers of N arose from infrared singularties 
appearing when internal fermion lines were restricted to momenta on the Fermi surface, 
similar to the Fermi surface enhancements discovered by S.-S. Lee for the problem of a 
Fermi surface coupled to a U(l) gauge field. These enhancements distinguish the present 
problem from that considered in Refs. 18 and 19: the Ising-nematic transition in a rf-wave 
superconductor. Formally, the latter problem is described by a field theory similar to that 
of the present paper: fermions with linear dispersion coupled via a Yukawa interaction to a 
scalar field </>. Also, in both problems we find a logarithmic divergence of velocity ratios in 
the infrared at order 1/A^ for the RG flows. However, for the d-wave superconductor, with 
Dirac fermions whose energy vanishes only at isolated "hot spots" , the 1/A^ expansion was 
found to be stable at higher loops. In contrast, for the present SDW problem, the fermion 
hot spots are connected to "cold" Fermi lines, and singularities associated with these lines 
lead to a breakdown in the naive 1/A^ counting. Because of this breakdown, the nature of 
the A^ — )■ oo limit of Eq. (2.1) remains unclear. 

Next, we examined the instability of the SDW metal to the onset of superconductivity 
near the quantum critical point in Section V. We found a strong tendency towards spin- 
singlet pairing, with pairing amplitude having opposite signs across a pair of hot spots. For 
the cuprate Fermi surface in Fig. 1 this includes (1^2 _y2 pairing, while for the pnictide Fermi 
surfaces this includes s_|__ pairing. This pairing instability was manifested in a log-squared 
divergence of the renormalization of the pairing vertex, arising from an interplay of the 
infrared singularities associated with the Fermi surfaces and the hot spot. This log-squared 
singularity cannot be resolved by the present RG approach, and other methods are needed 
to determine its consequences. An important problem for future research is to understand 
the feedback of the pairing fluctuations on the non-Fermi liquid singularities at the metallic 
hot spot. Clearly, superconductivity appears near the quantum critical point as T — !■ 0. The 
interesting question is the behavior above T^., involving the interplay between the metallic 
quantum criticality and the pairing fluctuations. 

In our discussion of the critical theory for the SDW transition in Section II, we noted 
that the field theory had emergent pseudospin SU(2) symmetries (Eq. (2.6)) containing the 
particle-hole transformation; note that the pseudospin rotations can be carried out indepen- 
dently on different pairs of hot spots. Given the strong instability towards d-wave pairing 
near the SDW critical point described in Section V, it is natural to examine the action of the 
SU(2) pseudospin symmetries on the (i-wave pairing order parameter. This was described 
in Section VI, where we found a similar log-squared enhancement of the susceptibility to a 
modulated valence bond solid (VBS) order parameter illustrated in Figs. 22 and 23. Notice 
that at short scales this ordering has an Ising-nematic character: this corresponds to the 
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breaking of a 90 degree rotation symmetry of the square lattice by the values of the bond 
order parameter in Eq. (6.10). It would be interesting if future work supports a connection 
between the ordering instability of Section VI, and the bond and Ising-nematic ordering ob- 
served in experiments^^"'^^. While the present analysis has focused exclusively on the vicinity 
of the hot spots, it is quite possible that strong coupling physics away from the hot spot 
could lock in a preference for commensurate values, such as Q^ = 0, in Eq. (6.10), leading 
to global Ising-nematic order. Also, it would be interesting to study the changes in the VBS 
ordering for the case of a SDW transition at an incommensurate ordering wavevector, like 
that found in the hole-doped cuprates. 

Finally, we note an interesting possibility for future theoretical work. Given the break- 
down of the 1/A^ expansion for the theory in Eq. (2.1) for the SDW critical point in a 
two-dimensional metal, other systematic methods of analyzing this field theory are clearly 
needed. Following Ref. 23, one possibility is to modify the (V0)^ term in Eq. (2.1) to k^^^cfP', 
where k is the momentum carried by </>. Then at the RPA level, we obtain a theory with 
z = 1 + X, and an expansion in small x appears possible. 
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Appendix A: RG computations 

In this appendix we give the details of our calculations in Sections II and III. 

1. RPA polarization 

We begin with the RPA polarization bubble, 

/rll rl^T 
^{G[{1 + q)Gi{l) + Gi{l + q)G[{l)) (Al) 

The two terms in brackets come from the two graphs in Fig. 4 with different directions of the 
particle flow. As discussed in Section II such graphs are equal by the emergent particle-hole 
symmetry. Thus, focusing on the contribution from ^ = 1, 

J [211)^ (^iriilr + Qr) -vi-{l + q)){ir]lr -V2-1) 

43 



We change variables to li = Vi ■ {I + q) , I2 = V2 ' ^ , and take the hmit r/ — )■ using the relation, 

1 P 

= irisgaMSix) (A3) 

X + zrj X 

which yields, 

U'=\q) = J^f -^ ( ^ + 7rzsgn(/, + qr)S{h)] (f + 7iisgn{lr)S{k)] + {q ^ -q) 

(A4) 
Evaluating the integrals over /i, I2, 

N f 

U^-\q) = / dlrSgn{lr + qr)sgn{lr) + (g ^ -q) (A5) 

QTTVxVy J 

Here, we've taken the principal value integral to be zero, as it would be if we used a particle- 
hole symmetric regularization. Otherwise, one can check that any terms generated by the 
pv integral are of the form iqr and are cancelled by the {q — )■ —q) term of Eq. A5. Now, 
subtracting the value of the polarization bubble at g = 0, we obtain, 

U'=\q) - U'=\q = 0) = -^^^ [dW (sgn(/. + g.)sgn(/.) - 1) + {q ^ -q) = ;^^^|g.| 

Siiv^Vy J 27rv^Vy 

(A6) 
which, taking into account contributions from the other hot spots, gives, 

Nv 
U{q) = U{q = 0) + - \q^\ (A7) 

2TTV^Vy 



2. Fermion self energy 

We next proceed to the self-energy of fermion ipl^^, Fig. 5, 

dLd'^l 



f dl d^I 

^i,.Ap) = r:,r;^' J -^G2{p - l)D{l) 



~ N ""' J (27r)3 ,^(p^ _ Q _ ^2 ■ (p - f) 7|/,| +P 

We take the limit 77 — )■ and use Eq. (A3). Moreover, we change variables, so that l± = V2-l 
and /|| is the momentum component along the Fermi surface of ip2 (i-e. perpendicular to -02). 
Then, 

V. / N 3 f dlrdl_i_dl\\ f P ■ n NT/, -- ^\ 1 

Si p = — — / : -_, + TTZsgn /^ - Vrmi. - «2 ■ p) 

N\v\ ] (27r)^ \l±-V2-p J 



j\ir\ + il + q 

(A9) 
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Thus, the imaginary part of S is given by, 

ImEi(p) = ^ / %gn{lr - pr)^==l==^ (AlO) 

N\V\ J StT \/l\lr\ + \V2-P\^ 

where we have performed the integral over l±, /y. Since, ImS(pT- = 0) = 0, 
ImSi(p) = -— / -^(sgn(/, -p,) - sgn(Z,))^== 



3 

-sg: 



n{pr) \^\/i\Pt\ + iv2-py - \v2 ■ P\) (All) 



2nN\v\'y 
On the other hand, the real part of S is given by, 

T3 V r ^ 3^2 -P f dlrdl\\ 1 1 



2^H-/ (2^)'y/^MT7j^l^-l + ^S + (^2-pT 



Changing variables to m = \/jIt + ^ 



II 



ReSl(p) = -;g^ / rf/|, / rf«-,2, ,^.,__ ., 



3(f;2-p1 / ci/|| .,,,-l/^^ (A13) 



27r'^N'y\v\ J |'02 ■ p| V 1^1 

The integral over /y is ultra-violet divergent. Cutting off the integral at |/||| = A, we obtain 
to logarithmic accuracy, 

ReSi(p) = -4|t^ log r^ (A14) 

7r^A'|w|7 |f2-pl 

Combining eqs. (All), (A14) we obtain the self-energy (3.26). 



3. Boson-fermion vertex 

Proceeding to the first correction in 1/A^ to the boson-fermion vertex. Fig. 6 

dlrd'' 



HI rt"^! 
5TlAp,q) = (rVV).., / -^G2{l+p)G,{l+p + q)D{l) (A15) 



Evaluating the matrix product. 



STM = ~j 



(2vr)3 ^2 ■ {1+ p) - I7]{lr + Pr) Vl ■ {I + P + q) - ir]{lr +Pr + Qr) 1\It\ + T^ 

(A16) 
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The integral (A16) is logarithmically divergent in the UV. To extract this divergence, we 
may set all external momenta to zero: 

A^ J (27r)3 {-vJx + Vyly - irilr){yJx + Vyly - ir]lr) 7|/r| +11 + 11 

The poles in ly coming from the two fermion propagators in Eq. (A17) are in the same half- 
plane; we may choose to close the ly integration contour in the opposite half-plane, picking 
up the pole from the bosonic propagator: 



UV 



6T{p,q) = - 



1 r dLdh 1 



''T""'a; 



N J (27r)2 {-vj^ - Wysgn{lr)^/-f\lr\ +ll){vJx - iVySgia{lr)^/-f\lr\ + 11) 

X — , ^ (A18) 



Changing variables to n = \/^y\l^'+l^, 



UV ^ r dh r dfu 

A^7loo2vr4| 2Ti'vlll + vlu^ 



^nP.<l) = — l ^1 --In--rr2 (A19) 



We now go to polar coordinates, Vxlx + iVyU = \v\pe 



-1 , 



UV J- I (^P 



5r(p,g) = — — -/ ^ de (A20) 

The integral over p is logarithmically divergent in the UV\ cutting off the integral at p ~ A, 

5r(p,g)'"= ^tan-i-logA (A21) 

TinJS a 



4. Boson self energy 

We now proceed to the 1/A^ corrections to the boson self-energy. Fig. 7. We first analyze 
the contribution of diagrams a),b) and c), which we label 5Ili. Utilizing the expression 
(2.15) for the fermion induced quartic coupling, we obtain, 

72^ 



<^nf(g) = ^|^r^-(g,-g,/,-/)Z^(/) 



dl d^l 

-^U^'^'^q. -q. I. -I) + r^'iq. i. -/, -q) + r'\q. i. -q, -i))d{i) 

(A22) 
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The first two terms in Eq. (A22) vanish (these terms correspond to the diagrams in Fig. 7 
a),b) ). Thus, only the diagram in Fig. 7 c) contributes, 



5n,(gr, q) = \qr\A{qr, q) + B{q,, q) (A23) 



with 



M,r, ^l-~^j i0Oi(l - ,)G'S + ,m) (A24) 

B(qr,q}^-^Y. f tr^\'r\G{(l - g)Gi{l + q)D(l) (A26) 

The quantity Alq^-, q) is logarithmically divergent in the UV . The coefficient of the di- 
vergence may be extracted by setting the external momenta and r to zero. Then, from 
Eq. (A16), we recognize, 

Mqr. q) '^ ^^ y 5V{p, q) = — tan"^ - log A (A26) 



T^VxVy ^ n-K a 



Now, let us evaluate B. We temporarily keep only the contribution from the hot spot pair 



,£=1/ -K J- / dl^dH 1 



B'=\qr.^ 



TTV^Vy J (27r)3 {vJx + Vyly " w 1 " ^ " z?7(/^ - q^)) 



X z l_d (A27) 

{-Vxlx + Vyly + V2- q- ivilr + qr)) {i\It\ +ll + E + r)' 



Note that the region |/t-| < |g^| does not contain any UV divergences. Thus, to compute the 
UV divergent part, we can confine our attention to the region |/t-| > \qr\. In this case, the 
two poles in ly coming from the fermion propagators in Eq. (A27) lie in the same half-plane; 
we may choose to close the ly integration contour in the opposite half-plane, picking up the 
pole from the bosonic propagator: 

,£=1/ -. UV i f dlrdlx 1 



B^-\q.,q-)^^^r 

1 1/ I 

(A28) 



T^v^Vy J\i^\>\q^\ (27r)2 vj^ _ ^y^sgn(/^) v/^T^TTfTT^ - Vi ■ q 

1 |/.| 

-Vxlx - iVySgn{lr) VtKtI +lx + r + V2- q2 ^/'y\lr\ +11 + r 



47 



Note that we may extend the integration over /^ in Eq. (A28) back to the whole real line 
without influencing the UV part of the result. Thus, 



B'=\qr,q) 



UV 



TTVxVy JO 



X 



dL 



1 



dlr 

27r y_oo 2vr (vj^ - ivy^^-flr + II + r - Vi ■ q) 
1 l^ 



{vj^ + iVy^J-ilQ + ll + r - V2 ■ q) 1^^lr + l\ 



c.c. (A29) 



It is convenient to change variables to u = v 7|/t| + ^^ + '^ 



B'=\qr.q) 



UV 



1 



TTV^VyY 



dL- 



du 



u'^ — ll — r 



27r Jy/l2^ 27r {vj^ - iVyU - vi ■ q){vJx + iVyU ~ V2 ■ q) 



+ c.c. 



(A30) 



The r in the lower limit of the integral over u may be dropped without influencing the UV 
behaviour. We now go to polar coordinates, v^lx + iVyU = |iT|pe*^, 



B'=\qr,q} 



UV 



'K{2'KVxVy-fy VxVy 



pdp 



'^-tan-la p2(^isin2 6)_Q,cos2 6') - ^ 



de 



tan ^ a 



(pe*^ — V2 ■ q){pe '^ — vi ■ q) 



+ c.c. 



(A31) 



The integral over p is quadratically divergent. Expanding the divergent part in g and r 

2 I^P 



B'=\qr.q) 



UV 



■TT— tan a 

pdp I dO 



irn VxVy J Jtan-i a 



a 



sin^ 9 — a cos^ 6^ 1 I 1 H — (f'l + V2) ■ qcos 9 



P 



+ \{{ih ■ q){v2 ■ q) + {{vi ■ q)^ + ih ■ q)^) cos29)j - 



VxVy r 



(A32) 



As usual, the term constant in q corresponds to a shift in the position of the critical point 
and will be dropped below. The term linear in g* vanishes under 9 ^ 71 — 9, i.e. Ix ^ —Ix 
(more rigorously, this term must vanish by symmetry, once the contributions from all 4 pairs 
of hot spots are summed). Finally, the term quadratic in g* and the term linear in r give 
logarithmic divergences. Cutting off the integral over p at p ~ A, 



5'='(g.,g1 = ^logA 



Txn^ 



tan 



a^ 



a 



a I + a^ 



a^ql I tan ^ 



a 



a 1 + a" 



r tan 



a 
(A33) 



Now, summing over the four pairs of hot spots, we restore rotational invariance. 



B{qr,q) 



im 



1 2 \ -1 1 

a + \ —z + a tan — 

a V a"' / a 



4 1 

g^ log Ah tan^^-rlogA (A34) 

Tin a 



We now compute the diagram in Fig. 7 d), which we label Sliu. This diagram is present 
already in the Hertz- Millis theory and, being momentum independent, leads only to a renor- 
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malization of r, 

dL(f!' 1 5ur f (fl 1 



.n.w^5./^z)(0S-4„./ 



2-k'^N 



(27r)3 (^|/^| + f2)2 7riV7 J (27r)2 p 

5 

■ur log A (A35) 



Now combining Eqs. (A23), (A26), (A34), (A35) we obtain the UV part of the correction 
to the boson propagator, Eq. (3.34). 



Appendix B: Violatations of large- A^ counting 

1. Boson- fermion vertex correction at three loops 

In this section we compute the vertex correction in Fig. 9. As shown in section IV, an 
attempt to evaluate this graph directly with bare fermion propagators results in infra-red 
divergences. To cure this problem, we dress the fermion propagators by the one-loop self- 
energy (3.26). For simplicity, we include only the imaginary part of the self-energy respon- 
sible for the dynamics. The frequency independent real part responsible for the logarithmic 
running of the velocity v will be ignored here. Thus, we use. 



-ijf9{^,vrk) + vf-k 



where 1 = 2, 2 = 1 and 



g(u, k) = sgn{uj){^-f\uj\ + P - \k\), Cf = —-— (B2) 

Then, the diagram in Fig. 9 is given by, 

^r,,.,t = -28Ar|^|ll^|i^G-(A:)G-(fc-/0Gr^(^-/2)G-(^^)G^(/0G}(/,) 

xD{k)D{k)D{k-k) (B3) 

The external fermions are taken to have hot spot index i = 1, while the fermions in the 
loop are taken to have i' = —1. As discussed in section IV, the contributions from i' = 2 
and i' = 4 are not enhanced in A^, while i' = 1 contributes a UV finite term of 0(1) when 
the external fermion momenta are chosen to lie on the Fermi surface. As we are mainly 
interested in corrections to mean-field scaling, we only retain UV divergent contributions 
below. Hence, all the external momenta of the diagram have been set to 0. Substituting the 
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one-loop corrected propagators (Bl), we obtain, 

ST = -28N f ^-^^-^^-^ ^ X i 

^^'^' J (27r)3 (27r)3 (27r)3 -t^^g^k^, v,.k)-v,.k'' -i^giK, O2 ■k)-v,-k 



-i^g{K - hr, vi- {k- h)) -V2- {k- k) -i^g{k^ - kr, h ■ {k - k)) - Vi ■ {k - h) 
^ ^ -D{h)D{h)D{h-h) (B4) 



-i^gihr, vi ■ h) +V2-li -i^gikr, ^2 • h) +vi-l2 

We may divide the spatial momenta into two groups: -Oi ■ k, V2 ■ k, V2 ■ h, Vi ■ h and vi ■ h, 
V2 ■ h- The singular manifold of the diagram is given by setting the momenta in the first 
group to zero and can be parameterized by the two variables in the second group. We begin 
by integrating over the first set of variables, picking up the contribution from the poles of the 
fermion propagators. As this integration is saturated at momenta of 0(l/iV), we can neglect 
the dependence of the boson propagators and fermion self-energies on these momenta. We 
then obtain the result in terms of an integral over the singular manifold. 



Due to the symmetry, G{1) = —G{—1), the contributions to the integral from kr > and 
kr < are equal. Now, changing momentum variables to vi ■ p, V2 ■ P, and integrating over 

V2 ■ k, Vi ■ I2, 



6r 



<t>il>2i'{ 



-7N- 



\v\ 



(v^n 



x'^yl JO 

dli 



dkr f d{vi ■ k)d{v2 ■ k)d{v\ ■ li)d{v2 ■ h) 

~2^ 



00 /"On 



27r 



00 /.O 

fcr 



dl 



2t 



2-K 



{2nY 
D{h)D{l2)D{h-l2) 



vi-l2=i>2-li=0 

1 



X 



-i^igihrA'i ■ h) - g{kr - hr, vi-{k- h))) + V2-k i%g{kr, vi-k) + V2 ■ k 



X 



-ijj{g{l2T,V2 ■ h) - gikr - kr, V2 ■ {k - k))) + Vi ■ k ij^g{kr, V2-k)+vi-k 

(B5) 

Now, performing the integral over Vi ■ k, V2 ■ k, 

■12 ''°° dkr r dhr r dl2r f d{Vl ■ h)d{v2 ' k) 



^^4>^2^l ^^ (,, „ ^3^2 



{vxVy)^cj Jo 2tt Jkr 27r Jk^ 2ti 

1 



(2vr) 



X 



g{kr, 0) + g{kT) + g{kT - kr, Vi ■ k) g{kr, 0) + gikr) + gikr - kr, h ■ k) 
xD{k)D{k)D{k-k) 



vi-l2=V2-h=0 
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Changing variables to /i,2t = ^ra^i,2, h,2y = VlKyi,2, 



^r,,,,t = -X(a) y^ ^ = X(«)logA, (B6) 



with 



2 ^oo 



X{a) = - i — + a] I dxi I dx2 I dyi I dy2- 



ISn^n \a J J^ J^ 7„^ 7„^ ^Xi + yf + ^Xi - 1 + yf - 2|y|i + 1 



1 1 1 

X r— i r^^^ X 



V^2 + yi + \/x2 -l + yl- 2\y\2 + 1 a;i + 4(a+")^^? 2;2 + 4(- + a)2i/2 

1 



2. Quartic vertex 

In this section we evaluate the five loop correcton to the boson four-point function shown 
in Fig. 18. We recall that by the particle-hole symmetry of our theory, diagrams with a 
reversed direction of the two fermion loops have the same value. We focus only on the 
diagrams where the fermions in the two loops come from opposite hot spots as these give 
a result, which is of 0{N^) and logarithmically divergent. To identify the coefficient of the 
logarithmic divergence we may set all the external momenta to zero. Then by rotational 
invariance each hot spot pair gives the same contribution. Moreover, we can also consider 
the diagram as in Fig. 18 but with fermions 1 and 2 interchanged. By refiection symmetry, 
this has the same UV divergence. Finally, we should be able to absorb the UV divergence 

^2 

into the coefficient of the quartic vertex 0^ , which specifies the spin structure, 

r-paia2a3a4 UV jj^ / ^0102^0304 , ^0103^0204 , raia4 ra2a3N r-p3333 {V,'R\ 

O 

and 

5r3333 = _4 . 6 . 2 . n . ^ . iV^ I '^^^^^^^^^^^^ - h)D{l2 - h) 

X G\{p^)G\{p^fG\{p^ - h)Glipi - /2)G}(j9i - h) 

X G^\p2)G2\p2?G^\p2 - k)G2\p2 - k)G-\p2 " /3) (B9) 

with 

S = tr{T\^T''T''T^T'^)tr{T^T^T^T^T^T'^) = 84 (BIO) 

We will used the same strategy for evaluating the integral (B9) as for computing the 
vertex correction in section Bl. The singular manifold in the present case is specified by 
vanishing pi, P2, Vi ■ li, V2 ■ h, Vi ■ I3 and can be parameterized by the three momenta V2 ■ h, 
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(B7) 



Vi ■ hi V2 ■ h- We will integrate explicitly over the first set of momenta and leave the result 
as an integral over the later three momenta. 

Let us call I{pit,P2t) the result of integrating over all momenta and frequencies in Eq. 
(B9), except pir and P2t- Then, using the particle-hole symmetry, G{p) = —G{—p), and 
the inversion symmetry, G~^{pr,p) = G^{pr, —p), we obtain I{pit,P2t) = H~Pit, —Pir) and 

I{PIt,P2t) =HP2r,Plr)- ThuS, 



^p3333 



2^u . 3^ . 7 . jY^ 



2 \ 5 poo 



dPlr f^^^ dP2T f dliT-dl2rdlsr 



2v^Vy I L O'TT I 9-rr I (9-rr\3 

d{vi ■ pi)d{v2 ■ Pi)d{vi ■ P2)d{v2 ■ P2)d{vi ■ h)d{v2 ■ li)d{vi ■ l2)d{v2 ■ k)d{vi ■ h)d{v2 ■ h) 



(2vr) 
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1 



-i^diPir, 0) + tTi . Pi (_i^^(p^^, 0) + iTa ■ pi)2 -i^g{p^^ - h^, V2 ■ h) + v^ ■ {pi - h) 



X 



'i^giPlT - kr, Vl ■ k) + V2 ■ ipi - h) -ijjg{.PlT - kr, V2 " k) + Vi ■ {pi - k) 



X 



-iw9{P2T, 0) - Vl ■ P2 {-i^9{p2r, 0) - V2 ■ P2)2 -^^^(p2r - kr, V2 " h) " ^^1 " (P2 " h 



1 



-ii!9{P2r - kr, Vl ■ k) - V2 " (p2 " k) -i^9{P2T " kr, ^2 ' k) - Vl " (P2 " k) 

xD{h)D{k)D{h-k)D{k-k) 
Integrating over Vi ■ k, V2 ■ k, Vi ■ k, 

5rf33 = -z2io.32.7-iV2 



(Bii: 



■17 poo 



dpir T'" dp2r 



{2v^Vy)^ ./n 27r ./_„-, 27r 



PlT 



00 /•P2 



PlT 



dkr 

~2^ 



00 j'P2 



PlT 



dl 



2t 



2-K 



00 l'P2 



PlT 



dl 



3r 



27r 



d{vi ■ Pi)d{v2 ■ Pi)d{vi ■ P2)d{v2 ■ P2)d{v2 ■ k)d{vi ■ k)d{v2 ■ k) 



{2ny' 



1 



X 



-'i'^giPiT, 0) + t/i ■ pi {-i^g{pir, 0) + tT2 ■ Pif 

1 1 



X 



^jj9{P2t, 0) - iTi ■ P2 {-i^g{p2r, 0) - iTa ■ pa)^ 



X 


-i^igipir - 


-kT,V2-k) + 9iP2r- 

1 


-kT,V2-k)) + vi 


ipl- 


-P2) 




-ij!^9{Pir - 


-kr^Vl -k) +9{P2r- 
1 


-kr^Vi -k)) +V2 


{Pl- 


-P2) 



-ij^iaiPlr - kr, V2 ■ k) + 9iP2r " kr, ^2 • k)) + ^1 " (Pl " P2) 

X D(/i)D(/3)i)(/i - k)D{k - /3)l,,.r.=...r.=...r3=o 
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(B12) 



Now, integrating over vi ■ pi, V2 ■ pi, 

|r?|5 rca 



6T 



3333 



i2'^ . 3^ ■ 7 ■ iV' 



dpir T'" dp2r f°° dkr 



{2VxVyJ JQ _„ ^„pj^ -.. jpj^^ 



5 '- 27r 
d({)i ■ P2)d{v2 ■ P2)d{v2 ■ h)d{vi ■ l2)d{v2 ■ h) 



27r /„, 271 
1 



(2vr) 



X 



1 



i49{P2r, 0) - Wi ■ P2 {-ijj9{P2r: 0) - ^2 " ^2)^ 



1 



'-i^{g{pir, 0) + 5'(/2r - Plr, ^^1 ■ ^2) + gikr " P2r, ^'l " ^2)) + ^2 " P2) 



X 



diij- I (JiI'St 



pi 



r 27r ip^^ 27r i2L(^g{j)^^^ q) + 5f(/i^ - pi^, V2 ■ k) + fi-l^ir - P2t, ^2 ■ /i)) - Vi ■ p2 

1 

i%{9iPlT, 0) + 5f(/3^ - Pi^, '02 • I3) + fl'(/3r - P2r, ^2 " /3)) - Vi ■ P2 



+ 



(iZiT 



P2t 



rf/ 



3t 



p^^ 27r ;_^ 27r _.£/ 



i^iaiPlT, 0) + 5'(/i^ - Pi^, V2 ■ h) + gihr - P2r, V2 " ^l)) + Vi ■ P2 

1 

i^igihr-PlT,V2 -h) +g{llT-P2T,V2-ll) +giPiT-kT,v2-h) +g{p2T-kT,v2-k)) 

^ dkr r dkr 1 



P2 



27r 



Plr 



27r _,•£/ 



iw{g{PlT, 0) + 5f(/3^ - Pir, V2 ■ k) + gikr - P2t, h ■ k)) +V1-P2 



i^{g{Plr - kT,V2 ■ k) + g{P2T -kT,V2 -k) + gikr -P1T,V2 ■ k) + g{kr -P2T,V2 ' k)) 



D{k)D{k)D{li - k)D{k - I. 



3/ \ vi-li=V2-l2=vi-h=0 



(B13) 



Observe that under k "^ ^3 the first term in the square brackets is invariant, while the 
second and third terms map into each other. Utihzing this fact and integrating over vi ■ p2, 

h ■ P2, 

u-ri3 



r"p3333 r)12 q2 



2'^ -3^ -T- N' 



dpir T'" dp2r P dkr r dkr 



{2vxVyCff Jo 2t\ Jq 2t\ Jp^^ 2 



in 



Pl-r 



27r 



/ 



d{v2 ■ k)d{vi ■ k)d{v2 ■ k) 



(27r) 



giPlr, 0) + g{p2r, 0) + g{kr " Plr, '^2 ' ^l) + gikr " P2t, h " k) 
1 



X 



{g{PlT, 0) + g{p2r, 0) + g{kT - Plr, Vi ■ k) + gikr " P2r, C'l " ^2))^ 

''- dkr 1 

Plr 2vr g{p^^^ 0) + g{p2r, 0) + g{kr - Plr, V2 ■ k) + gikr " P2r , ^2 " ^3) 

"'^dkr^ 1 

-00 27r g(^l^^ - p^^^ £,2 . k) + g{kr - P2r, V2 " ^l) + g{Plr " kr, V2 ' k) + g{P2r " kr, V2 " ^3 

D{k)Dik)Dik - k)D{k - /3)|,,,-;=...r.=.,r3=o (B14) 

We now introduce dimensionless variables, p2r = xpir, kr = ViPir, i'2 ■ k = ^IPirZi-, 



+ 
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— * — * 

vi-k = y/7PirZ2, ^ ' k = y/lPirZs- Then, 

^p3333 ^ lN'Y{ah r^ = N'Y {ah log A (B15) 

2 Jo Pit 

with 

rr? / -I \ /*! /"OO /»00 /"OO /"OO /"OO 

^<"' = -27^U + "J io'^'/ **"/ '"" L"" J-J'' L"'' 



1 



X 



X 



;i + V^ + Vz/2 - 1 + 4 + ^1/2 - X + 2| - 21221)=^ 

dy-i — — 

1 1 + ^ + a/i/3 -1 + zl + a/z/s - X + 2| - 2|23| 

+ / dy'i 



X 



a/i/i - 1 + 2^ + V'yi - a; + ^f + Vl - 1/3 + ^1 + Va; - ^3 + z| - 2|2;i| - 2|23| 
1 1 1 



"" lz/2 - 1/31 + \C- + «)2(^2 + zl) + l(a2 - J,);,2;^3 ^^^^) 



3. Pairing vertex 

This appendix will describe the direct evaluation of the pairing vertex correction in 
Eq. (5.5). We first attempt to perform the calculation using bare fermion propagators, 



5V 



Vip'itp'i 



— 3/i f dlrdl±dl\\ 



N\v\^J (27r)3 7|/.| + /i + /jJ l^-y,.k,-^^il^-k 

1 



It) 

V\ 



X 



\v\ 



where we've introduced variables l± = V2 ■ I, l\\ = ^ij{v2)ilj- For simplicity, let us choose 
kir = k-ir = u > 0. We now perform the integral over l±. For |/t-| < u both poles in the 
fermion propagators are in the same half-plane and we can pick up just the pole from the 
bosonic propagator. In the opposite regime, |/t-| > u, we get contributions from both the 
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bosonic and fermionic poles. Thus, 



3^ 

N\v\' 






2r]J^ 2n J 27r \l^ - ^lll . [k, + k_,) -fl^ + q + {v2 ■ h^ 



1 1 



(B19) 



The contribution from the bosonic pole in Eq. (B17) gives an expected logarithmic diver- 
gence, 

S'°'Ty^,^, ~ -—f-r. log — ^ (B20) 



y^tv^t ~ ^ 1^1^ log - — ^ 



On the other hand, the contribution from the fermionic poles in Eqs. (B18),(B19) gives 
a much stronger infra-red singularity. If we set the total momentum of the fermion pair 
ki + k_i to zero, then 

47rA'r7|f2 ■ /ci| \\V2 ■ ki 
with _ 



/(a) = / dx-^== (B22) 

If the total pair momentum is non- vanishing, in particular, if -|'i?2"(^i + ^-i)| ^ ('i'^2'^i)^;7i^, 
then. 



As usual, we cure the strong infra-red divergences by using a one-loop dressed fermion 
propagator (Bl). Then, 



i\\v\ 



X 



V2- (1 + k_i) + i—f-g{lr + A;-ir, ^^i ■ (/ -F fc.i)) 



For simplicity, we take the external fermion momenta to lie at the hot spots, ki = A;_i = 0. 
Moreover, as before, we choose the external frequencies, kir = /c_i,- = a; > 0. Switching to 
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variables, l±, /y, we perform the integral over l±. As we saw above, the contribution from 
the pole in the bosonic propagator could be calculated without dressing the fermion Green's 
function and was of 0{1/N) - we drop this piece below. On the other hand, as we will see 
the contribution from the poles in fermionic propagators is of (9(1) in A^. Moreover, since 
/_L ~ 0{1/N) at these poles, we may ignore the dependence of the fermion self-energy on 
/_L, which gives, -Oi ■ / = ^Jty^h- In this manner, we obtain. 



6^ f^ dlr f dill 1 



Cf\v\ J^ 27r J 27r 7/,- + In 
1 



(B24) 



We now perform the integral over /^. This integral is convergent in the ultra-violet. However, 
when w — )■ 0, it is logarithmically divergent in the infra-red. This infra-red divergence comes 
from the region 7/^ ^ /jj. Changing variables to 7/,- = xl?,, we obtain, 

3/i f°° dill f dx 



^'l\v\cf Jo h J| X + 1 ^a: + (^)2-| + ^^ + ( J^)^ + f - ^ 

(B25) 
For I?, ^ '-fu, performing the integral over x to logarithmic accuracy. 



6^a f°° dl\i -, ( ^\\ \ l^ot 2 [ ^ 




,2 



<^ry^t^t ^ —Tj^TT-^—r. / ^ log — = - r 2^n l^S — (^26) 

7r^7|w|c/(a^ + 1) i^^ /|j \-iUj ) 7r(a^ + 1) \^u ; 

4. Density vertex 

In this appendix, we compute the one-loop renormalization of the density-wave vertex, 
shown in Fig. 20b), 

f dH 

(jro^^t(A;i, fc„i) = 3^y ^D{l)G\{h - l)G2\k^i - I). (B27) 

If we ignore the effects of Fermi-surface curvature, G{1) = —G{—1), and Eq. (B27) reduces to 
its counterpart in the superconducting channel with k_,i -^ — /c_i. In the present calculation, 
we will keep the effects of the Fermi-surface curvature using a propagator, 

GUI) = — ^ ^ ^ (B28) 

-^gil^,v^.l) + vf.l + inlyir 

Here, we ignore any dressing of the curvature by the interactions. 

For simplicity, we set external momenta to zero and choose ki^ = —k_ij. = a; > 0. As in 
Appendix B3, we introduce variables l± = V2- 1, l\\ = ^ij{v2)ilj- Proceeding as in Section V, 
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we keep only the contribution to the integral (B27) from the Fermi liquid regime, 7/^ ^ l?i- 
Then, 

_ 3/i /• dill f dlr f dh 1 

1 1 

Performing the integral over l_j_, 

.p 3/i fdlll f^V^dlrZilll) K ,^... 

^1^1 i 27r i^ 27r /j / ^(/||);f 

^r + I 2m 

Notice that the Fermi-surface curvature is present in the denominator of Eq. (B30). This 
is in contrast to the corresponding calculation in the superconducting channel, where the 
Fermi- surface curvature drops out. Performing the integral over l^, 

ffo**, = -^ /" ^^ '°S 'j (B31) 

where we have ignored terms subleading in /y in the numerator of the logarithm. Recall, 
Z(Z||) ~ A^|i7|/||. Hence, for /y <^ {mu/N\v\) ' the /^ integral is cut-off in the infrared by 
the external frequency and the Fermi surface curvature may be neglected. On the other 
hand, for /y ^ {muj/N\v\) ' the integral is cut-off by the curvature. By comparison, in the 
superconducting channel the integral is cut-off by the external frequency in both regimes 
resulting in a stronger enhancement. Notice that the cross-over scale {muj/N\v\) ' is much 
larger than the infra-red cut-off of the /y integral ^JjJJ. Evaluating the integral over /y to 
leading logarithmic accuracy, 

^ro^^t = - , ,7^,, log^ f — ) (B32) 

^^ 37r(a^ -I- 1) \1^ J 
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